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ABSTRACT 


This  thesis  developed  a  propagation  model  that  can  be  used  to  investigate  the  perfor¬ 
mance  of  coherent  and  incoherent  beam  configurations  subjected  to  diffraction  and  thermal 
blooming  effects.  The  propagation  model  also  provides  an  additional  damage  assessment 
parameter — power-in-the -bucket — as  well  as  the  irradiance  in  the  target  plane  for  any  laser 
wavelength.  The  propagation  model  results  have  target  irradiance  patterns  similar  to  those 
obtained  from  WaveTrain,  although  the  peak  irradiance  differs  by  ~30%  due  to  differences 
in  thermal  blooming  modeling.  In  addition,  the  propagation  model  can  optimize  certain 
beamlet  configurations’  peak  irradiance  using  a  genetic  algorithm.  Based  on  all  the  sim¬ 
ulation  results,  a  13-beamlet  configuration  is  deemed  to  have  highest  peak  irradiance  and 
power-in-the -bucket  than  all  the  tested  multiple  beam  configurations. 
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CHAPTER  1: 
Introduction 


Over  the  past  several  decades,  the  United  States  has  researched  and  developed  high-energy 
lasers  for  different  platforms  [1],  Although  Directed  Energy  (DE)  weapons  are  complex  to 
design  and  costly,  they  deliver  their  payload  at  the  speed  of  light  [1],  However,  the  dwell 
time  of  a  DE  weapon  imposes  an  engagement  time  constraint  that  reduces  the  DE  weapon’s 
capability  of  handing  swarm  and  high-speed  targets.  Therefore,  a  short  dwell  time  is  desired 
for  DE  weapons  during  target  engagement.  An  output  of  ~  1  MW  is  required  to  kill  many 
kinds  of  target  in  less  than  eight  seconds  based  on  an  estimate  developed  in  Chapter  6. 

The  entire  thesis  is  be  based  on  Solid  State  Lasers  (SSL)  due  to  their  deep  magazine 
capability,  relatively  low  cheap  cost  per  kill,  technological  maturity  and  other  factors  to  be 
discussed  in  Chapter  2.  The  current-state-of-the-art  for  a  single  solid-state  laser  source  is 
only  tens  of  kilowatts  [2].  One  approach  to  build  a  ~  1  MW  output  power  laser  system  is 
to  spatially  combine  arrays  of  available  high-power  fiber  lasers.  In  order  to  evaluate  the 
performance  of  spatially  combined  laser  beams,  atmospheric  modeling  is  essential  as  the 
laser  beam  propagates  from  the  DE  weapon  through  the  atmosphere  to  the  target. 

Due  to  time  constraints,  the  implementation  of  this  propagation  model  only  includes  diffrac¬ 
tion  and  thermal  blooming  effects  for  various  laser  beam  configurations.  In  Chapter  2,  the 
thesis  starts  with  an  overview  of  DE  weapons  and  end  with  the  different  types  of  available 
laser  sources.  In  Chapter  3,  the  thesis  discussed  various  techniques  to  combine  laser  beams. 
In  Chapter  4,  the  thesis  discussed  the  atmospheric  effects  that  a  propagating  laser  beam 
would  encounter.  In  Chapter  5,  the  thesis  modeling  tools  is  discussed.  In  Chapter  6,  the 
thesis  validated  the  model  with  analytical  computations  and  simulations  result  from  another 
propagation  modeling  tool,  WaveTrain.  In  Chapter  7,  the  thesis  discussed  and  concluded  the 
performance  differences  between  single  source  megawatt  output  lasers  and  those  obtained 
using  beam  combination  techniques  during  target  engagement. 
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CHAPTER  2: 

Overview  of  Directed  Energy  Weapons 


Unlike  kinetic  energy  projectile  weapons,  DE  weapons  use  high  irradiance  electromagnetic 
waves  of  a  certain  wavelength  to  inflict  damage  onto  the  intended  target.  The  main  damage 
mechanism  from  lasers  is  the  ability  to  deposit  huge  thermal  energy  onto  the  target.  The 
thermal  energy  deposited  can  result  in  melting,  penetrating,  and/or  malfunctioning  of  the 
target’s  hardware.  The  actual  degree  of  damage  depends  highly  on  the  laser’s  power,  dwell 
time,  and  beam  quality.  Table  A.l  in  the  Appendix  provides  an  estimate  of  the  power  level 
required  to  counter  certain  targets. 


2.1  History  of  DE 

In  1954,  the  Office  of  Naval  Research  funded  Charles  Townes  to  develop  the  laser.  Charles 
Townes  envisioned  that  lasers  would  be  used  as  oscillators  for  optical  communication. 
Gordon  Gould  proposed  to  the  U.S.  Advanced  Research  Projects  Agency  (ARPA)  that  a 
laser  could  be  used  to  measure  an  enemy  target’s  distance  or  act  as  a  designator  for  missiles. 
He  also  concluded  that  a  laser  beam  could  be  10000  times  brighter  than  the  sun,  an  intensity 
that  might  "trigger  chemical  reactions  or  perhaps  ultimately  cause  nuclear  fusion"  in  nuclear 
missiles  [3].  The  first  laser  was  developed  by  Theodore  Maiman  from  Hughes  Research 
Laboratories  in  1960  [3].  The  United  States  Air  Force  (USAF)  awarded  Maiman  and 
various  military  labs  contracts  to  develop  lasers  for  missile  guidance  and  communications. 
The  first  optical  laser  was  based  on  a  ruby  gain  medium  and  hence  the  output  power  was 
limited  to  tens  of  watts  due  to  the  ruby’s  material  properties  and  lack  of  good  thermal 
management.  The  Department  of  Defense  (DOD)  realized  the  potential  of  the  laser  as  an 
operational  military  system  from  these  inventions  and  invested  heavily  in  its  research  and 
development.  The  Air  Force  Chief  of  Staff,  General  Curtis  LeMay,  commented  on  March 
28,  1962,  that  "directed  energy  weapons  would  be  able  to  transmit  energy  across  space  with 
the  speed  of  light  and  bring  about  the  technological  disarmament  of  nuclear  weapons."  [3]. 
By  mid- 1965,  C.  Kumar  N.  Patel  at  Bell  Labs  developed  a  200  W  continuous  wave  gas  laser. 
Subsequently,  Hughes  Research  Laboratories  developed  a  1.5  kW  with  a  10  m  oscillator 
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and  a  54  m  amplifier. 


The  key  breakthrough  for  high-energy  lasers  came  from  addressing  the  heat  transfer  problem. 
Arthur  Kantrowitz  and  Ed  Gerry  at  the  Avco  Everett  Research  Laboratory  near  Boston 
addressed  the  laser’s  thermal  issue  by  using  a  rocket-like  approach:  burning  fuel  to  generate 
hot  CO2  and  mixing  with  nitrogen  before  expanding  through  special  nozzles  at  supersonic 
velocity  [3].  The  highly  excited  CO2  molecules  created  the  necessary  population  inversion 
conditions  for  lasing.  In  the  1970s,  the  USAF  Airborne  Laser  laboratories  developed  the  first 
aerial  platform  laser  system,  based  on  a  400  kW  CO2  laser.  From  the  1980s  to  the  present, 
many  laser  systems  such  as  the  Mid-infrared  Advance  Chemical  Laser  (MIRACL),  Airborne 
Laser  (ABL),  Laser  Weapon  System  (LaWS),  etc.,  have  been  successfully  developed  and 
tested.  After  the  market  appearance  of  the  first  commercial  fiber  laser  in  late  1980s,  higher 
power  single-mode  fiber  lasers  have  been  continuously  developed  for  machining.  These 
high-power  fiber  lasers  are  potential  sources  for  DE  weapons  due  to  their  high  wall  plug 
efficiency.  Finally,  after  decades  of  research  and  development  on  high-energy  military 
lasers,  DE  weapons  with  the  capability  of  countering  certain  targets  at  greater  distances  are 
ready  to  be  installed  in  Navy  ships. 


2.2  Advantages  of  DE  Weapons 

The  main  advantage  of  DE  weapons  over  kinetic  weapons  is  their  engagement  speed  and 
precision.  Kinetic  weapons  require  time  for  their  projectiles  to  reach  the  target,  and  hence, 
engagement  would  be  at  closer  range  when  compared  to  DE  weapons,  which  engage  at  the 
speed  of  light.  This  fast  engagement  speed  also  gives  DE  an  edge  when  handling  rapidly 
maneuvering  targets.  In  general,  kinetic  energy  weapons  cause  collateral  damage,  unlike 
DE  weapons,  which  allow  the  user  to  precisely  select  the  aimpoint.  The  cost  of  target 
engagement  for  DE  weapons  is  lower  than  that  of  tactical  missiles  as  only  electrical  energy 
is  required  for  laser  operations.  Although  the  initial  setup  costs  for  DE  are  high,  the  cost  per 
shot  and  the  other  advantages  mentioned  make  DE  weapons  a  worthy  long-term  investment. 
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2.3  Disadvantages  of  DE  Weapons 


DE  weapons  require  line-of-sight,  and  thus  obscured  targets  cannot  be  engaged.  Targets 
that  are  temporally  obscured  will  increase  the  engagement  time  requirements.  Due  to  these 
engagement  time  requirements,  DE  can  only  handle  a  limited  amount  of  swarm  targets, 
unlike  kinetic  energy  weapons.  Atmospheric  conditions  can  also  adversely  affect  DE 
weapons’  range  when  compared  to  kinetic  energy  weapons. 


2.4  DE  Source  Technologies 

The  source  for  DE  weapons  can  be  gas,  chemical,  solid  state,  or  free  electron  lasers.  Some 
of  these  laser  sources  have  been  successfully  deployed  by  the  United  States.  Gas  technology 
was  briefly  covered  in  the  history  section.  The  following  sections  will  focus  on  chemical, 
solid  state,  and  free  electron  lasers. 

2.4.1  Chemical  Laser 

Chemical  lasers  rely  on  chemical  reactions  for  their  pumping  energy  and  more  specifically, 
typically  operate  on  molecular  transitions,  except  in  the  case  of  atomic  iodine.  Chemical 
lasers  produce  near-to-mid  infrared  spectrum  wavelengths.  In  these  devices,  mixed  chemical 
vapors  flow  through  the  gain  region  in  a  direction  that  is  transverse  to  the  laser  beam 
axis.  Due  to  chemical  reactions,  toxic  and  combustible  chemicals  are  produced,  and  hence 
residual  gases  must  either  be  neutralized  or  collected  for  later  disposal.  The  logistic  handling 
of  chemicals  must  also  be  considered  in  operational  scenarios.  To  maximize  the  laser’s 
efficiency,  a  special  nozzle  is  used  to  control  the  flow  conditions  within  the  gain  region. 
Although  dealing  with  chemical  lasers  can  be  cumbersome,  chemical  lasers  have  been 
successfully  tested  for  military  applications.  The  U.S.  Navy  developed  the  first  megawatt 
class  laser  known  as  the  MIRACL,  which  was  in  operation  from  1980  to  2000.  MIRACL 
used  the  chemical  deuterium  fluoride  to  produce  a  4  /mi  laser  and  was  integrated  with  the 
SeaLite  Beam  director  for  operational  testing.  Due  to  MIRACL’s  wavelength,  it  was  prone 
to  absorption  during  atmospheric  propagation.  In  the  1990s,  the  USAF  developed  another 
megawatt-class  laser  known  as  the  ABL.  The  ABL  used  six  chemical  oxygen  iodine  1.3  /mi 
lasers  within  a  modified  Boeing  747  airplane.  In  2010  field  testing,  the  ABL  successfully 
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intercepted  and  destroyed  ballistic  missiles  during  their  initial  boost  phase  [1].  However, 
the  program  was  discontinued  due  to  economic  and  practical  issues. 

2.4.2  Solid-State  Laser 

Certain  materials  that  consist  of  crystals  or  glass  with  specific  atomic  ions  imbedded  are 
excellent  for  the  laser  gain  medium.  The  specific  atomic  ions  must  have  a  long  radiative 
decay  lifetime  and  favorable  absorption  band,  allowing  pumping  light  to  be  effectively 
absorbed.  SSLs  use  optical  pumping  to  create  population  inversion  conditions  for  lasing. 
The  optical  pumping  can  be  provided  by  another  laser,  a  semiconductor  diode,  or  a  flashlamp. 
There  are  various  types  of  SSL  structures  but  only  fiber  and  slab  SSLs  have  been  used  as 
DE  weapons. 

2.4.2.1  Slab  Solid-State  Lasers 

In  a  slab  SSL,  the  gain  medium  is  a  large  plate-like  structure.  The  DOD’s  Joint  High 
Power  SSL  (JHPSSL)  program  developed  a  DE  prototype  known  as  the  Maritime  Laser 
Demonstration  (MLD)  that  leveraged  slab  SSL  technology.  In  2010,  the  MLD  project  group 
coherently  combined  seven  15  kW  slab  SSLs  to  create  a  1.064  /rm  laser  with  an  output 
power  of  105  kW  [1],  During  field  testing,  the  MLD  successfully  tracked  and  engaged 
small  boats  in  a  marine  environment.  In  May  2011,  a  Northrop  representative  stated  that 
they  "could  build  the  first  unit  of  a  full-power  engineering  and  manufacturing  development 
(EMD)  version  of  the  MLD  within  four  years,  if  the  Navy  could  find  the  resources  to  fund 
the  effort"  [1],  However,  given  the  current-state-of-the-art,  slab  SSLs  cannot  be  scaled 
to  megawatt  power  levels  as  the  slab  is  not  readily  cooled  due  to  low  surface  to  volume 
ratio  [1], 

2.4.2.2  Fiber  SSL 

In  a  fiber  SSL,  the  gain  medium  is  a  flexible  fiber  cable  that  is  more  readily  cooled  than 
a  slab,  due  to  an  improved  surface  to  volume  ratio.  A  fiber  SSL  consist  of  a  laser  diode 
pump,  a  wavelength-selective  coupler  that  allows  the  pump  wavelength  to  enter  the  fiber 
gain  medium  without  disturbing  the  signal,  an  amplifier  spliced  into  the  fiber  gain  medium. 
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and  an  optical  isolator.  High  power  fiber  lasers  are  widely  used  in  material  cutting  and  for 
DE  weapons.  One  of  the  most  prominent  full  operational  fiber  lasers  system  is  the  Laser 
Weapon  Systems  (LaWS),  which  has  been  deployed  on  the  USS  Ponce  since  2014. 

2.4.23  LaWS 

The  LaWS  project  started  in  2010  [4]  and  was  awarded  to  Kratos  Defense  and  Security 
solutions  by  the  Naval  Surface  Warfare  Center.  LaWS  was  designed  to  be  integrated  with  the 
existing  Phalanx  Close-in-Weapon  System  (CIWS)  with  the  intent  to  disable/jam  Electro- 
Optical  (EO)  sensor  suites  and  to  counter  UAVs/  EO  guided  bombs  with  the  assistance  of 
radar.  LaWS  combined  six  5 .5kW  IPG  Photonics  fiber  lasers  to  produce  a  total  output  power 
of  33  kW  with  a  beam  quality  of  17  [1],  The  wavelength  of  each  fiber  laser  is  1.06  //m  and 
is  near  an  atmospheric  transmission  window  (1.045  /mi).  To  ensure  LaWS  effectiveness, 
it  was  tested  with  different  intended  targets  under  various  combat  scenarios.  In  June  2009, 
LaWS  accurately  engaged  five  UAVs  at  China  Lake  in  southern  California  [1],  In  May  2010, 
LaWS  accurately  engaged  four  UAVs  at  1  nautical  mile  away  from  San  Nicholas  Island  [1], 
During  the  same  tests,  LaWS  also  demonstrated  the  ability  to  destroy  rigid-hull  inflatable 
boat  materials  at  about  half  a  nautical  mile  and  also  disabled  EO/infrared  sensors.  Between 
July  and  September  2012,  LaWS  accurately  engaged  three  UAVs  when  the  system  was  on 
an  Arleigh  Burke  class  destroyer.  LaWS  was  installed  on  the  USS  Ponce  for  a  one  year  trial 
as  of  August  2014  and  was  declared  as  an  operational  asset  in  September  2014  [1],  If  all 
the  testing  goes  well,  the  U.S.  Navy  could  deploy  laser  weapons  operationally  in  the  near 
future. 


2.4.3  Free  Electron  Laser 

A  Lree  Electron  Laser  (EEL)  consists  of  free  electrons  in  vacuum  as  the  gain  medium,  a 
radio-frequency  linear  accelerator  (linac)  as  the  pumping  source,  alternating  magnetic  field 
as  an  undulator,  laser  resonator  mirrors,  and  an  output  coupler  [5].  Unlike  conventional 
lasers  that  rely  on  population  inversion  in  a  gas  or  solid-state  material,  a  EEL  uses  relativistic 
electrons  from  a  linac,  which  pass  through  an  undulator,  causing  them  to  oscillate  and 
therefore  produce  light.  During  these  oscillations,  the  electrons  radiate  in  a  relativistic  dipole 
radiation  pattern  that  maximizes  in  a  direction  perpendicular  to  the  oscillation  direction.  The 
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electrons’  relativistic  speed  causes  the  oscillation  frequency  to  shift  from  a  lower  frequency 
in  their  rest  frame  to  a  higher  frequency  in  the  laboratory  frame.  The  radiation  frequency 
depends  on  the  electrons’  kinetic  energy  and  the  period  of  the  alternating  magnetic  field. 
Mirrors  are  placed  at  opposite  ends  of  the  laser  cavity  to  reflect  a  portion  of  the  radiated 
energy  back  through  the  alternating  magnetic  fields.  The  reflected  light  stimulates  additional 
radiation  that  produces  a  strong  optical  beam  within  and  beyond  the  laser  cavity.  Existing 
FELs  have  wavelength  ranges  from  0.06  nm  to  3  mm  [6]  with  scalable  output  power  up  to 
1  GW  for  pulsed  lasing  and  10  W  for  continuous  lasing.  The  current-state-of-the-art  for 
high  average  power  FELs  is  14  kW  at  Jefferson  National  lab.  The  FEL’s  wavelength  can 
be  tuned  for  optimal  atmospheric  transmission  and  has  excellent  beam  quality  suitable  for 
DE  applications.  FEL  technology  is  thought  to  be  scalable  to  MW  class  since  there  is  no 
conventional  substrate  to  damage.  The  energy  of  spent  electrons  can  be  recovered  to  increase 
wall  plug  efficiency.  Currently,  FELs  tend  to  be  large(~20  m),  heavy,  and  expensive.  During 
FEL  operation,  additional  radiation  shielding  is  required  due  to  Bremsstrahlung  radiation 
produced  by  the  electrons.  Around  2010,  the  U.S.  Navy  had  a  plan  known  as  Innovative 
Naval  Prototype  (INP)  to  develop  a  ~  1 00 kW  class  FEL,  but  the  INP  was  shelved  in  ~  2012 
due  to  the  immaturity  of  FEL  technology,  and  the  technology  advancement  in  fiber  laser, 
the  U.S.  Navy  may  eventually  reassess  FEL  technology  to  achieve  ~MW  power  level  single 
laser  source  system. 
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CHAPTER  3: 

Overview  of  Beam  Combination  Techniques 


A  weapon’s  effectiveness  on  the  battlefield  can  be  measured  using  the  probability  of  damage. 
Fragmentation  and  blast  warheads  increase  the  probability  of  damage  by  ensuring  sufficient 
coverage  area.  For  DE  weapons,  the  probability  of  damage  depends  on  beam  power,  beam 
quality,  dwell  time,  slant  range,  weather,  atmospheric  conditions,  etc.  The  beam  power 
refers  to  the  total  amount  of  power  contained  within  the  laser  beam  as  it  exits  the  DE 
weapon.  The  beam  quality  parameter  characterizes  its  focusability,  which  affects  the  laser’s 
irradiance  and  spot  size  on  the  target’s  surface.  Dwell  time  refers  to  the  time  required  for 
the  laser  to  damage  the  target. 

For  DE  weapons  to  be  considered  lethal  for  many  applications,  they  need  to  destroy  targets 
within  seconds  of  engagement  and  thus  they  require  a  laser  beam  power  from  hundreds  of 
kilowatts  to  megawatts  with  good  beam  quality;  the  actual  output  power  needed  depends 
on  the  target  and  the  atmospheric  conditions,  and  the  laser  system  itself.  The  continuous 
advancement  of  laser  technologies  might  eventually  lead  to  the  development  of  a  single  laser 
source  that  has  hundreds  of  kilowatts  or  even  megawatts  of  output  power.  Current  solid-state 
lasers  are  largely  limited  by  thermo-optic  or  other  non-linear  effects  that  set  a  limit  on  the 
output  power.  Although  fiber  lasers  have  large  surface  to  volume  cooling  area,  they  are 
still  vulnerable  to  thermally  induced  damage  above  a  threshold  power.  A  single  mode  fiber 
laser’s  output  power  will  saturate  at  several  kilowatts  unless  there  is  a  major  breakthrough 
in  the  fiber’s  mode  control  and  better  thermal  solutions  [2].  An  alternative  approach  to 
building  high-power  laser  systems  is  to  combine  arrays  of  commercially  available  high 
power  fibers.  Over  the  past  several  decades,  researchers  have  advanced  the  state  of  laser 
beam  combination  technology.  Currently,  it  is  possible  to  design  a  laser  system  with  tens 
of  kilowatts  of  output  power  and  good  beam  quality  using  an  array  of  laser  beams. 

Beam  combining  techniques  for  laser  arrays  can  be  characterized  in  three  broad  classes: 
incoherent,  coherent,  and  spectral  beam  combining.  These  combination  techniques  can 
have  either  a  tiled- aperture  or  a  filled- aperture  architecture.  The  tiled- aperture  architecture 
has  several  laser  beams  aligned  with  optics,  allowing  them  to  be  bundled  at  the  output 
aperture  of  a  DE  weapon,  whereas  the  filled-aperture  architecture  has  these  several  beams 
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overlaid  to  cover  the  entire  beam  director  mirror,  as  illustrated  in  Figure  3.0.1.  To  compare 
beam  combination  techniques,  one  could  consider  the  laser’s  irradiance  on  the  target  or  the 
power  in  the  bucket. 


Figure  3.0.1:  An  example  of  tiled-aperture  architecture  (left)  with  individual 
laser  beamlets  represented  by  different  colors  and  filled-aperture  architec¬ 
ture  (right)  showing  a  single  output  beam. 


3.1  Incoherent  Beam  Combination 

The  working  principle  of  the  incoherent  beam  combining  (IBC)  technique  is  to  focus  all  the 
beamlets  onto  the  target  without  accounting  for  their  individual  phases.  The  IBC  technique 
can  be  implemented  using  either  a  filled-aperture  or  a  titled-aperture  architecture.  For  the 
filled-aperture  architecture,  the  laser  beamlets  are  combined  in  a  fiber  that  degrades  the 
beam  quality  significantly  and  thus  is  not  suitable  for  DE  weapons.  For  the  tiled- aperture 
architecture  shown  on  the  left  in  Figure  3.0.1,  the  laser  beamlets  are  bundled  side-by- 
side.  The  beamlets  may  or  may  not  have  the  same  wavelength;  even  those  with  the  same 
wavelength  may  have  random  wavefront  phase  or  polarization  relative  to  one  other.  The 
beamlets  must  be  aligned  so  that  they  will  overlap  at  the  same  spot  on  the  target.  The  total 
irradiance  produced  by  IBC  will  be  N  times  the  individual  beamlet  irradiance  where  N  is  the 
number  of  beamlets,  and  the  overall  combined  beam  waist  diameter  is  inversely  proportional 
to  each  beamlet’ s  diameter.  Figure  3.1.1  shows  the  intensity  pattern  of  six  beamlets  at  the 
beam  director’s  output  aperture.  Figure  3.1.2  shows  the  far-held  pattern  of  a  beam  formed  by 
the  IBC  technique  after  the  six  beamlets  propagate  through  the  atmosphere  [7] .  Figure  3.0. 1 
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and  Figure  3.1.1  do  not  have  any  central  beamlet  as  the  beam  director’s  secondary  mirror 
occupies  that  position  as  shown  on  the  right  in  Figure  3.1.4. 


Figure  3.1.1:  Intensity  of  six  beamlets  at  the  output  aperture  [7],  The  X 
and  Y  axes,  which  represent  the  beamlets'  position  at  the  output  aperture, 
are  arbitrary  units.  The  Z  axis,  which  represents  the  intensity,  is  normalized 
to  the  beamlet's  peak  intensity  at  the  output  aperture. 
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Figure  3.1.2:  Far-field  pattern  of  a  beam  formed  by  IBC  [7],  The  X  and 
Y  axes,  which  represent  the  combined  beam’s  position  at  far-field,  are  arbi¬ 
trary  units.  The  Z  axis,  which  represents  the  intensity,  is  normalized  to  the 
beamlet’s  peak  intensity  at  the  output  aperture. 


The  IBC  technique  has  numerous  advantages  over  coherent  and  spectral  combining.  This 
approach  does  not  require  the  individual  beamlets  to  have  narrow  line  widths  or  to  have  the 
same  phase  or  polarization.  Without  considering  the  atmospheric  turbulence  effects,  the 
spot  size  of  a  Gaussian  laser  beam  propagating  in  free  space  is  given  as 

",b)  =  "4  +  (i) )  ■  <3'° 

where  Zr  =  nw02/  A  is  the  Rayleigh  length  and  wa  is  the  spot  size  at  the  focus,  as  shown  in 
Figure  3.1.3. 
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radius  w(z). 


To  obtain  a  small  spot  size  and  thus  high  irradiance  on  the  target,  each  beamlet  must 
have  good  optical  beam  quality,  and  its  diameter  must  be  sufficiently  large  at  the  output 
aperture.  According  to  Equation  3.1,  a  small  wQ  spot  will  result  in  shorter  Rayleigh  length. 
Therefore,  the  system  requires  multiple  large  beam  expanders  or  telescopes.  In  theory,  IBC 
laser  systems  can  be  scalable  to  any  power  level  required  for  DE  applications,  but  the  beam 
director  radius  is  a  function  of  each  beamlet  diameter(a),  R  Beam  Du-  ~  Nl/2a,  as  shown  in 
Figure  3.1.4  [8]. 


Figure  3.1.4:  Beamlets  bundled  side  by 
is  shown  on  the  left  and  the  actual  bear 


side  at  the  beam  director  aperture 
i  director  is  shown  on  the  right. 
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3.2  Coherent  Beam  Combination 


The  working  principle  of  the  coherent  beam  combining  (CBC)  technique  is  to  get  all  the 
beamlet’s  individual  phases  to  be  the  same  by  using  phase  correction  mechanisms.  The 
total  irradiance  produced  by  CBC  will  be  N2  times  the  individual  beam  irradiance,  and 
the  overall  combined  beam  waist  diameter  is  approximately  inversely  proportional  to  the 
overall  diameter  of  the  array.  Figure  3.2.1  shows  the  far-held  pattern  of  a  beam  formed  by 
the  CBC  technique  from  six  beamlets  at  the  beam  director’s  output  aperture  as  shown  in 
Figure  3.1.1.  Figure3.2.1  also  shows  side  lobes  that  are  formed  by  constructive  interference 
of  the  individual  beamlets  [7]. 


Figure  3.2.1:  Beam  formed  by  CBC  technique  that  uses  tiled-aperture  con¬ 
figuration  [7],  The  X  and  Y  axes,  which  represent  the  combined  beam’s 
position  at  far-field,  are  arbitrary  units.  The  Z  axis,  which  represents  the  in¬ 
tensity,  is  normalized  to  the  beamlet's  peak  intensity  at  the  output  aperture. 


The  CBC  technique  is  subdivided  into  active  and  passive  techniques. 
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3.2.1  Active  Coherent  Beam  Combining 

Active  coherent  beam  combining  techniques  can  have  either  a  filled-aperture  or  tiled- 
aperture  architecture.  For  the  tiled- aperture  architecture,  a  master  oscillator  power  amplifier 
(MOPA)  and  electronic  phase  feedback  system  are  used  to  alter  the  individual  laser  beamlets’ 
phases  as  shown  in  Figure  3.2.2. 


Master 

s. V 

Oscillator 

'  ^  Phase 
Modulators 


Expander 


Output 


Phase 

Control 


Figure  3.2.2:  Active  CBC  configuration  using  MOPA  and  electronic  feed¬ 
back. 


To  provide  a  coherent  input  to  each  individual  laser’s  amplifier,  the  MOPA  output  must  be 
polarized  with  a  narrow  line  width  spectrum.  The  electronic  phase  feedback  system  can 
actively  adjust  each  individual  beamlet’s  phase  to  compensate  for  any  phase  differences 
between  the  beamlets. 

The  phase  feedback  system  can  also  be  used  as  a  small  angle  beam  steerer  to  assist  aiming 
by  varying  the  beamlet  phases  to  change  the  constructive  interference  pattern  on  the  target. 
This  change  in  aiming  would  be  due  to  changes  in  the  combined  beam’s  wavefront.  This 
beam  steering  method  is  similar  to  the  working  principles  of  active  phased  array  antennas 
in  radar. 

The  fiber  length  differences  between  the  beamlets  are  very  crucial  for  the  computation  of 
the  phase  feedback,  which  is  the  input  for  the  path  length  adjustors.  Since  the  phase  of 
beamlets  can  be  adjusted,  it  is  possible  to  compensate  for  low  atmospheric  turbulence. 

Any  active  CBC  that  uses  a  tiled- aperture  configuration  would  introduce  side  lobes  that 
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would  decrease  the  maximum  peak  irradiance  on  the  target  if  there  are  gaps  within  the 
beamlets,  as  shown  in  Figure  3.2.1.  In  addition,  variation  between  the  individual  amplitudes 
of  the  beamlets  could  increase  these  side  lobes. 


The  tiled-aperture  CBC  technique  cannot  achieve  more  than  75%  irradiance  when  compared 
with  a  perfect  Gaussian  beam  due  to  its  side  lobes  [8].  For  the  coherent  tiled  aperture 
beam  combiner,  the  larger  separation  between  the  beamlets  and  the  non-uniform  illuminated 
subapertures  reduce  the  amount  of  irradiance  on  the  target.  A  waveguide  combiner  can 
theoretically  be  used  to  compensate  for  the  irradiance  drop.  In  a  square  or  rectangular 
waveguide,  an  input  beam  reproduces  itself  after  travelling  for  a  distance  known  as  the 
Talbot  length  given  by 


(3.2) 


where  n  is  the  refractive  index  inside  the  waveguide,  and  h  is  the  internal  dimension  of  the 
waveguide.  If  the  waveguide  terminates  before  the  Talbot  length,  these  beams  are  separated 
into  multiple  beams  that  are  then  coherently  combined  into  a  single  output  beam.  If  a  larger 
separation  between  the  entry  beams  is  desired,  a  tapered  waveguide  must  be  designed.  Lock¬ 
heed  Martin  Corporation  has  been  conducting  waveguide  coherent  combining  experiments 
for  a  tiled- aperture  architecture  [8].  Northrop  Grumman  Corporation  researchers  have  also 
used  another  method,  involving  a  diffractive  optical  element  rather  than  a  waveguide,  to 
produce  high  beam  quality  beams  for  a  coherently  combined  tiled- aperture  architecture  [8]. 
Although  the  fill-factor  of  the  output  aperture  has  increased,  the  electronic  steering  and 
automatic  aberration  ability  of  the  beams  have  decreased. 


3.2.2  Passive  Coherent  Beam  Combining 

The  passive  CBC  techniques  use  optical  feedback  from  the  fibers  to  lock  all  the  beamlets’ 
phases.  Current  examples  are  the  fiber  ring  and  the  self-Fourier  cavity.  Unlike  the  active 
CBC  technique,  the  passive  CBC  techniques  do  not  require  narrow  line  width  output.  This 
allows  a  single  fiber  to  operate  at  a  higher  output  since  the  onset  of  damaging  nonlinear 
effects  is  delayed  for  broader  spectrum  beams.  However,  there  are  still  technical  issues 
associated  with  the  scalability  of  the  passive  CBC  techniques,  and  hence  this  thesis  will  not 
focus  on  them. 
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3.2.2.1  Fiber  Ring 


This  technique  uses  a  spatial  filter,  which  includes  a  focusing  lens  and  pinhole,  to  select 
those  longitudinal  modes  that  interfere  with  constructive  laser  combination.  Fiber  lasers 
have  many  longitudinal  modes  due  to  their  length  and  large  line  width,  and  hence  the 
probability  of  finding  the  appropriate  longitudinal  mode  for  constructive  interference  is 
higher.  The  fiber  ring  CBC  technique  can  be  mathematically  described  by  [8] 


K 

=I>< 

7=1 


,  Aj  J 


(3.3) 


where  K  is  the  number  of  fibers  used,  n  is  the  fiber  core  refractive  index,  Aj  is  the  mode 
amplitude,  and  Lj  is  the  length  of  jth  fiber.  The  feedback  signals  to  the  fiber  amplifiers  are 
proportional  to  \E2\  and  are  maximized  when  the  wavelengths  are 

nL; 

~T  =  mj  ,  (3.4) 

aj 

where  mj  is  an  integer.  Currently,  Northrop  Grumman  Information  Technology  Sector 
researchers  have  successfully  combined  16  fibers  using  this  technique  [8].  The  passive  ring 
CBC  technique  produces  good  laser  beam  quality  and  is  relatively  simple  as  it  does  not  re¬ 
quire  any  phase  locking  system  or  polarization-maintaining  fiber.  Since  it  can  automatically 
select  any  of  the  longitudinal  modes,  it  can  compensate  for  any  thermal  and  mechanical 
perturbation.  This  technique  enables  higher  output  power  scalability  since  it  can  operate 
with  broadband  fibers,  which  are  inexpensive.  However,  this  technique  requires  a  lens  that 
might  distort  the  wavefront  due  to  the  high  heat  load  imparted  from  the  high-power  sources. 


3.2.2.2  Self-Fourier  Cavity 

Unlike  the  Fiber  Ring  CBC  technique,  the  Self-Fourier  technique  relies  on  a  Fourier  lens  to 
spatially  separate  the  beamlets  before  reflecting  them  back  to  the  input  sources  for  coherent 
combination.  All  these  inputs  sources,  which  reside  in  the  input  plane  one  focal  length 
away  from  the  Fourier  lens,  are  pointing  along  the  optical  axis  and  hence  are  centered  on 
the  optic  axis  in  the  Fourier  plane,  which  is  one  focal  length  beyond  the  lens. 
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Consider  a  beamlet  located  at  a  distance  d\  away  from  the  optical  axis.  After  passing 
through  the  Fourier  lens,  it  is  collimated  with  a  propagation  angle  given  by 


d=f«n,  (3.5) 

where  F  is  the  focal  length  of  the  Fourier  lens.  Due  to  the  small  propagation  angle,  the 
phase  along  the  Fourier  plane  is 

a t b  =  kxsin(6 )  =  kx— ,  (3.6) 

F 

where  k  is  the  wavenumber,  and  x  is  the  distance  from  the  optical  axis.  If  x  is  a  multiple  of 

di  =  f,  (3.7) 

d\ 

the  phase  will  be  a  multiple  of  2n  and  will  constructively  interfere  in  the  Fourier  plane  with 
patterns  do  apart.  If  do  is  equal  to  d  \ .  then 


</i  =  VlF,  (3.8) 

and  the  interference  pattern  at  the  Fourier  plane  will  be  similar  to  the  original  pattern  of 
sources  at  the  input  plane. 

If  the  Fourier  plane  is  configured  to  coincide  with  the  input  plane,  which  can  be  implemented 
by  cutting  the  Fourier  lens  in  half  and  adding  a  partially  reflective  mirror,  the  feedback  pattern 
will  theoretically  be  coupled  back  to  the  input  sources  when  the  beams  are  coherent,  as 
shown  in  Figure  3.2.3.  As  the  beamlet’s  power  increases,  temperature  variation  across  the 
surface  of  the  lens  and  mirror  may  exist  as  cooling  can  only  occur  at  their  edges  to  prevent 
interfering  with  the  optical  path.  This  temperature  variation  will  affect  the  geometric 
characteristic  of  the  lens  and  thus  resulting  in  beam  distortion. 
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3.3  Spectral  Beam  Combination 

The  spectral  beam  combining  (SBC)  technique  uses  a  dispersive  optical  element,  such  as 
a  diffraction  grating  or  prism,  to  spatially  combine  beams  of  different  wavelength.  This 
technique  is  similar  to  the  wavelength-division  multiplexing  technique  used  in  optical 
communications.  The  SBC  technique  requires  each  beam’s  direction  and  wavelength  to  be 
controlled.  The  implementation  of  the  SBC  technique  is  easier  than  CBC  as  controlling 
wavelength  is  simpler  than  phase.  Laser  resonators  that  are  combined  using  the  SBC 
technique  share  the  same  output  mirror.  These  beams  are  parallel  and  are  focused  at 
the  desired  range.  For  a  particular  laser  element,  the  incident  angle  that  the  beam  hits 
the  dispersive  optical  element  determines  the  oscillating  frequency  in  the  case  where  the 
grating  element  is  part  of  the  cavity.  The  bandwidth  of  the  laser  gain  mediums  must  be  able 
to  support  all  these  oscillation  frequencies. 

Together  with  the  U.S.  Air  Force  Research  Laboratory,  MIT  Lincoln  Laboratories  developed 
a  SBC  architecture  that  uses  a  lens  to  focus  the  parallel  laser  inputs  [9].  Each  of  these  inputs 
is  separated  by  distance  Xj  and  then  focused  by  the  lens  on  the  dispersive  optical  element, 
as  shown  in  Figure  3.3.1. 
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Figure  3.3.1:  SBC  architecture  developed  by  MIT  Lincoln  Laboratories  and 
U.S.  Air  Force  Research  Laboratory.  Source:  [9]. 


The  diffraction  angle  caused  by  the  grating  material  depends  on  the  incoming  beam’s  wave¬ 
length.  An  incident  plane  wave  with  wavelength  A  propagating  in  the  direction  orthogonal 
to  the  grating  grooves  will  be  diffracted  according  to  [9], 


sin(jS)  -  sin(or) 


in  A 

~T 


(3.9) 


where  m  is  the  diffraction  order;  a  and  fi  are  the  incident  angle  and  reflected  angle  of  the 
wave,  respectively,  at  the  grating  surface;  and  /  is  the  grating  period.  Taking  the  derivative 
of  equation  3.9  with  respect  to  incident  angle  a  will  result  in 


Ad  _  /  cos(a) 

A  a  m 

The  fiber  separation  is  approximately 


(3.10) 


Xj  =  A  ax  f. 


(3.11) 
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where  f  is  the  focal  length.  Substituting  Atr  from  equation  3.10  into  equation  3.1 1, 


IX ;  cos(ff) 

AA  = - 7  / 

mf 


(3.12) 


From  equation  3.12,  the  spatial  separation  of  the  fibers  determines  the  difference  in  each 
fiber’s  operating  wavelength,  and  hence  the  wavelength  will  not  be  evenly  distributed  for 
constant  fiber  spacing  because  of  the  cos  (a)  term.  In  addition,  the  diffracted  beam  will 
have  an  angular  spread  that  is  larger  than  that  of  the  incident  beam  due  to  each  fiber’s  finite 
line  width.  The  divergence  angle  deviation  is  given  by 


A  A*  m 

l  COS(yS) 


(3.13) 


where  Ad*  is  the  bandwidth  of  an  individual  fiber.  Due  to  this  additional  angular  spread, 
the  beam  quality,  M 2,  becomes  worse  as  compared  to  a  single  fiber  laser  source. 


M 2 


AA*mna>0 
2/dcos(y0 ) 


2 


(3.14) 


where  oj0  is  the  beam  waist.  The  power  scaling  of  SBC  is  limited  by  the  grating’s  thermal 
capacity,  fiber  gain  bandwidth,  and  power  limits  for  fibers  with  narrow  linewidth.  To 
illustrate  the  beam  quality  obtainable  using  the  SBC  technique,  consider  a  source  comprised 
of  Gaussian  beams  with  wavelength  of  A  =  1  /mi.  The  source’s  beam  waist  is  assumed  to 
be  u>o  =1-5  cm  with  Ad*  =  1  x  10-9m  and  the  grating’s  blaze  spacing,  /,  is  assumed  to  be 
1  yum.  Using  these  assumptions,  M 2  =  23.6  according  to  equation  3.14  when  coSyS  =  1. 
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CHAPTER  4: 

Overview  of  Atmospheric  Propagation 


The  Earth’s  atmosphere  contains  many  different  constituents  that  vary  based  on  location  and 
time.  The  presence  of  these  constituents  introduces  atmospheric  absorption  and  scattering 
effects  that  reduce  a  DE  weapon’s  effectiveness  by  attenuating  the  optical  beam  energy. 
In  addition,  the  density  and  temperature  of  the  Earth’s  atmosphere  both  fluctuate  and  thus 
can  introduce  turbulence,  further  reducing  a  DE  weapon’s  effectiveness.  Because  a  DE 
weapon’s  laser  wavelength  is  typically  in  the  infrared  region,  this  chapter  will  only  focus  on 
the  atmospheric  effects  that  impact  these  wavelengths. 


4.1  Atmospheric  Absorption 

The  atmosphere’s  constituents  can  be  classified  into  two  main  categories:  molecules  and 
aerosols.  As  shown  in  Table  4.1,  the  Earth’s  atmosphere  is  mainly  comprised  of  nitrogen, 
oxygen,  argon,  carbon  dioxide,  and  water  vapor.  The  nitrogen  and  oxygen  molecules 
are  linear,  homonuclear  diatomic  molecules  that  have  symmetrical  bonds.  Due  to  these 
symmetrical  bonds,  their  dipole  moments  will  not  absorb  any  propagating  laser  beam  that 
is  in  the  infrared  spectrum.  Unlike  the  nitrogen  and  oxygen  molecules,  the  remaining 
molecules  have  asymmetrical  bonds  causing  their  dipole  moments  to  absorb  radiation  in 
certain  parts  of  the  infrared  spectrum.  Thus,  maximum  energy  can  only  be  deposited  on 
the  target  if  the  laser  beam’s  wavelengths  are  not  within  the  absorption  regions  of  these 
molecules.  The  aerosol  absorption  effects  on  propagating  laser  beams  are  similar  to  the 
molecular  effects,  but  these  effects  depend  substantially  on  the  aerosols’  size  and  content. 
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Content 

Composition 

Nitrogen 

78% 

Oxygen 

21% 

Argon 

0.9% 

Carbon  dioxide 

0.03% 

Water  Vapor 

0-4% 

Table  4.1:  Typical  earth’s  atmosphere  contents.  Adapted  from  [10]. 


4.2  Atmospheric  Scattering 

Additionally,  laser  beams  are  often  scattered  during  propagation  by  molecules  and  sus¬ 
pended  aerosols  in  the  Earth’s  atmosphere.  The  nature  of  these  scattering  effects  depends 
on  the  optical  beam’s  wavelength  and  the  size  of  the  scattering  agents.  For  wavelengths  that 
are  smaller  than  ~5  //m.  the  probability  of  scattering  an  optical  beam  is  inversely  propor¬ 
tional  to  the  fourth  power  of  the  optical  beam’s  wavelength.  This  type  of  scattering,  known 
as  Rayleigh  scattering,  occurs  when  the  scatterer  is  much  smaller  than  the  wavelength.  The 
nitrogen  and  oxygen  molecules  in  the  Earth’s  atmosphere  are  the  primary  causes  of  this 
Rayleigh  scattering.  For  example,  the  shorter  wavelengths  from  the  sun  are  scattered  by 
the  Earth’s  atmosphere,  which  is  why  the  sky  appears  blue  during  the  day.  The  deflection 
of  photons  tends  to  be  random  Rayleigh  scattering.  As  the  size  of  the  scattering  agent 
increases,  full  Mie  scattering  occurs.  Mie  scattering  is  more  directional  than  Rayleigh 
scattering  as  it  increases  in  the  forward  direction.  The  most  common  Mie  scattering  agents 
in  the  Earth’s  atmosphere  are  fog,  smoke,  and  dust  particles.  Mie  scattering  plays  an  impor¬ 
tant  role  in  laser  beam  propagation  since  these  beams  are  in  the  infrared  spectrum.  When 
the  scattering  agents’  size  becomes  ten  times  larger  than  the  optical  beam’s  wavelength, 
non-selective  scattering  occurs.  Mist  and  fog  appear  to  be  white  due  to  Mie  scattering 
effects.  Figure  4.2.1  shows  the  three  scattering  processes  [11]. 
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Nonselective  Scattering 


Figure  4.2.1:  Scattering  probabil  ity(%)  versus  wavelength. 

Adapted  from  [11], 


The  combined  effects  from  absorption  and  scattering  for  an  optical  beam  propagating 
through  a  linear  medium  can  be  modeled  by  Beer’s  law, 


P(z)  =  P0e~ez  ,  (4.1) 

where  P(z)  represents  the  power  of  the  optical  beam  at  z  distance  away  from  the  source  and 
P0  is  the  output  laser  power  from  the  DE  weapon.  The  variable  e  is  known  as  the  extinction 
coefficient,  given  by 


6  —  <ym  +  (Xa  “t  Pm  “t  Pa  >  (4.2) 

where  a  and  P  represent  absorption  and  scattering  coefficients,  respectively,  and  subscripts 
a  and  m  represent  molecular  and  aerosols  contributions,  respectively.  Figure  4.2.2,  4.2.3, 
4.2.4, 4.3.1  illustrate  the  absorption  and  scattering  coefficient  due  to  molecular  and  aerosols. 
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Figure  4.2.2:  Aerosol  absorption  extinction  coefficients  for  different  wave 
length  generated  using  MODTRAN. 


Figure  4.2.3:  Aerosol  scattering  extinction  coefficients  for  different  wave 
length  generated  using  MODTRAN. 


26 


Molecular  Absorption 


Figure  4.2.4:  Molecular  absorption  extinction  coefficient  for  different  wave 
length  generated  using  MODTRAN. 


Figure  4.2.5:  Molecular  scattering  extinction  coefficient  for  different  wave 
length  generated  using  MODTRAN. 
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4.3  Atmospheric  Turbulence 


Fluctuations  in  the  density  and  temperature  of  the  Earth’s  atmosphere  can  cause  the  air’s 
index  of  refraction  to  vary,  thus  changing  the  beam’s  propagation  pattern,  affecting  the 
direction  and  focus  of  the  beam.  The  index’s  variation  due  to  temperature,  wavelength,  and 
pressure  can  be  modeled  as  [12] 


p 

n  *  1  +  77.6  x  10“6[1  +  7.52  x  KT3/T2]-  ,  (4.3) 

where  n  is  the  refractive  index,  A  is  the  optical  beam  wavelength  in  pm.  P  is  the  pressure 
in  mbar,  and  T  is  absolute  temperature  in  Kelvins. 

Due  to  this  variation,  the  propagating  laser  beam  can  spread,  affecting  the  aim  point  as 
shown  in  Figure  4.3.1.  The  intended  laser  spot  on  the  target  in  Figure  4.3.1  is  in  the  center, 
but  the  final  spot  location  is  different  due  to  the  effect  of  turbulence. 
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Figure  4.3.1:  A  laser  beam  propagating  in  a  turbulent  atmosphere  in  X-Z 
plane  (left).  The  laser  spot  on  the  target  plane  (right).  Both  figures  are 
from  a  wave  propagation  simulation.  The  color  bar  axis  represents  relative 
irradiance. 


Turbulence  development  is  a  nonlinear  process  that  can  be  described  by  the  Navier-Stokes 
equations.  However,  the  solutions  to  the  Navier-Stokes  equation  are  unstable  for  finite  per¬ 
turbations  at  large  Reynolds  numbers.  Instead,  Kolmogorov’s  theory  based  on  a  statistical 
method  is  used  to  predict  turbulence  using  the  C2  value  [13].  The  C2  value  represents  the 
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amount  of  turbulence  in  the  atmosphere  that  the  beam  is  propagating  through.  Strong  tur¬ 
bulence  exhibits  a  C2  value  of  10-13m  ^  or  more,  while  a  typical  weak  turbulence  exhibits 
a  C“  value  of  10 or  less.  Due  to  less  influence  of  surface  heating,  the  C“  value  at 
higher  altitudes  tends  to  decrease. 

The  Fried  parameter,  r0,  measures  how  well  the  laser  beam  can  maintain  coherence  through 
the  atmosphere  that  has  varying  refractive  index.  For  a  spherical  wave,  the  ra  parameter  is 
given  by 


r0  *  0.33— ^  ,  (4.4) 

(/C„2)5 

where  /  is  the  distance  from  the  target.  When  turbulence  is  sufficiently  strong,  ra  is  used  to 
determine  the  target  spot  size(Equation  illustrated  in  Chapter  3)  instead  of  the  beam  director 
diameter. 


4.4  Thermal  Blooming 

As  the  high-power  laser  propagates  through  the  atmosphere,  some  parts  of  its  energy 
are  absorbed  by  molecules  and  aerosols.  This  absorption  process  causes  the  air  to  have 
temperature  gradients  that  ultimately  change  the  air’s  density.  The  heated  air  expands 
out  radially,  decreasing  the  air’s  density  and  hence  reducing  the  index  of  refraction.  This 
reduction  in  the  index  of  refraction,  modeled  by  Equation  4.3,  is  equivalent  to  placing  a 
diverging  lens  along  the  beam’s  propagation  path,  reducing  the  beam’s  peak  irradiance  on 
the  target.  This  effect  is  known  as  thermal  blooming. 

Before  deriving  the  thermal  blooming  equations,  the  beam  propagation  must  first  be  ana¬ 
lyzed  using  the  paraxial  wave  equation  [14], 

=  +  3TV1^  +  (■ ik6n >  (4-5) 

oz  2  2k 

where  e  is  the  extinction  coefficient  in  Equation  4.  2;  72*  V2  ifj  represents  the  beam’s  transverse 
diffraction;  and  (ikdn)tj/  describes  the  effects  of  the  refraction  index  variation,  cm,  caused 
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by  thermal  blooming  and  potentially  other  effects  like  turbulence.  The  general  solution  to 
Equation  4.5  is  given  by 


=  Aeikz  ,  (4.6) 

where  A  =  \A\e~u,>t  ,  k  =  2jtIa,  oj  is  the  angular  frequency,  and  z  is  the  distance  in  the 
direction  of  propagation.  The  refraction  index  variation,  6n,  caused  by  thermal  blooming 
is  given  by 


on 


Sn  Krr, 

— AT  =  — 
6T 


nn  -  1 


-AT 


(4.7) 


where  na  is  the  ambient  refractive  index,  T„  is  the  ambient  temperature,  and  AT  is  the 
temperature  change  in  that  particular  portion  of  the  atmosphere  due  to  absorption. 

This  temperature  change  due  to  thermal  blooming  is  described  using  the  isobaric  equa¬ 
tion  [14], 


dAT 

dt 


a 

pcP 


1  -  (v  •  V  JA T  + 


(4.8) 


The  a/PcpI  term  describes  the  heating  effect  from  the  laser,  where  a  is  the  absorption 
coefficient  in  Equation  4.2,  p  is  the  density  of  air,  Cp  is  the  constant-pressure  specific  heat 
capacity,  and  I  is  the  laser’s  irradiance.  The  (v  •  V  ±)AT  term  describes  the  convection 
effects  in  the  air,  where  v  is  the  wind  velocity,  and  Vj_  is  the  transverse  gradient  operator. 
The  (*7pcpV]_ )  AT  term  describes  the  conduction  effects  in  the  air,  where  K  is  the  thermal 
conductivity.  Equation  4.8  assumes  that  the  atmospheric  pressure  remains  constant  when 
the  propagating  high-power  laser  increases  the  temperature. 

In  a  situation  where  there  is  little  or  no  wind,  the  conduction  term  in  Equation  4.8  dominates. 
In  this  scenario,  the  thermal  blooming  effect  can  be  modeled  by  a  dimensionless  scaling 
factor,  e~Dc.  The  Dc  parameter  is  given  by  [15] 
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Dr  = 


njPaz' 


(4.9) 


2  nKn0a2 

where  az  <<  1,  nr  is  dn/dr  in  Equation  4.7,  P  is  the  laser  power,  and  a  is  the  beam  radius. 


In  windy  situations,  the  thermal  blooming  effect  depends  on  the  wind  speed  as  depicted  by 
the  second  term  on  the  right  in  Equation  4.8.  The  wind,  which  blows  across  the  propagation 
path,  causes  the  air’s  temperature  to  increase  downwind.  This  increase  in  temperature 
causes  the  air’s  refractive  index  to  decrease,  and  the  beam  bends  in  the  direction  of  the 
increasing  refractive  index  upwind.  When  the  wind  blows  in  the  transverse  x  direction  and 
dkr/dt  =  0(i.e.,  at  steady  state),  the  phase  shift  due  to  thermal  blooming  is  given  by  [14] 


r  ^  yi  \  CY  k  7  r  X 

<p(x,y,z)  =  ik  8ndz  =  — °— - Idx'  .  (4.10) 

Jo  To  pCpV  J — oo 

The  laser  irradiance  pattern,  for  a  given  (p ,  due  to  thermal  blooming  can  be  estimated  by 
solving  Equation  4.5  with  the  ray-optics  limit  using  perturbation  theory  [14].  The  solution 
to  Equation  4.5  can  be  approximated  to  the  first  order  by 


I(x,  y,  z)  =  4exp 


2  2 
X  x  +y 

2—e  °2 

a 


1  +  erf 


(4.11) 


where  IQ  is  the  irradiance  of  the  unperturbed  beam,  and  Nc  is  the  beam  irradiance  distortion 
parameter  accounting  for  the  wind  effects  given  by 


( na  -  1  )aPz2 
nn0pCpT  va3 


(4.12) 


The  larger  the  Nc  parameter,  the  more  thermal  blooming  deforms  the  propagating  beam’s 
irradiance  pattern  into  a  crescent  shape  shown  in  Figure  4.4.1. 
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Figure  4.4.1:  The  laser  irradiance  pattern  on  the  target  due  thermal  blooming 
with  Nc  =  1.  The  color  map  represents  the  irradiance  magnitude(i.e.,  dark 
red  represents  the  highest  intensity  while  light  blue  represents  the  lowest 
intensity). 
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CHAPTER  5: 
Simulation  Methods 


Another  mathematical  form  that  is  equivalent  to  Equation  4.5  is  the  time  invariant  Helmholtz 
equation  given  by 

[V2  +  k2n2]  E(r)=0,  (5.1) 

where  E(r)  is  the  complex  electric  field  of  the  propagating  laser  beam,  r  =  xi  +  yj  +  zk, 
k  =  2,1  / a  is  the  wavenumber,  and  n  is  the  index  of  refraction. 

5.1  Diffraction  Modeling 

A  solution  to  Equation  5.1  within  the  paraxial  approximation  is  given  by  the  Fresnel 
diffraction  integral  [16] 

JkAz  poo 

E(r2)  =  — -  /  E(ny  2E|r2-ri1  dn  ,  (5.2) 

iAAz  J- oo 

whereri  =  jcii+  r/ij  is  the  source  plane  vector,  r2  =  x2i+  i)2]  is  the  target  vector  coordinates, 
and  A z  is  the  propagating  distance  between  the  source  and  target  plane  as  illustrated  by 
Figure  5.1.1. 


y2 


Figure  5.1.1:  Optical  wave  propagation  coordinate  system.  Source:  [12], 
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It  is  very  tedious  to  solve  Equation  5.2  analytically  and  hence  it  would  be  advantageous  to 
use  a  computer  instead.  For  efficient  computer  simulation,  Equation  5.2  can  be  expressed 
in  Fourier  transform  terms  as  [12] 


gi*Az  fjfelnl  .  , 

E(r2)  = - — - Tmvx)e1^^-}  |fl  ,  (5.3) 

iAAz 

where  fi  =  (fx\  =  X2/a\z  ,  fy i  =  l,2/,i-\z)  is  the  source  spatial  frequency  with  grid  spacing 
of  5 1.  For  on-axis  phase,  the  e,kAz  term  can  be  ignored.  Equation  5.2  can  be  expressed 
as  another  Fourier  transform  equation,  which  performs  a  convolution  between  the  source 
plane  field  and  the  free  space  amplitude  spread  function,  as  [12] 


E(r2)  =  E(ri) 


aikAz 


iAAz 


-eVA1 


JkAz  ,1 

e  Ini2 1 


E(r2)  =  T~x  |r{E(n)}  x  T  j  — 


where  the  free  space  amplitude  spread  Fourier  transform  function  is 


(5.4) 


T  jei'MV'2fe|ri|2J  =  e^ze-inAAz\h?  _  (5.5) 

Substituting  Equation  5.5  into  Equation  5.4 

E(r2)  =  T~x  \T (E(ri)}  x  >  (5.6) 

Since  Equation  5.4  performs  two  Fourier  transforms,  Equation  5.3  is  desired  for  faster 
computational  speed.  Although  Equation  5.4  is  not  as  efficient,  it  compensates  by  intro¬ 
ducing  an  additional  parameter  which  provides  grid  spacing  flexibility  for  the  observation 
plane  [12].  The  equation  implemented  for  the  newly  developed  propagation  model  is  based 
on  Equation  5.4.  The  propagation  model  is  implemented  in  Matlab  to  simulate  incoher¬ 
ent  and  coherent  beam  combination  considering  diffraction,  thermal  blooming,  turbulence, 
absorption  and  scattering  effects. 

Since  the  computer  simulation  only  samples  E(r2)  at  a  fixed  interval,  it  is  difficult  to  obtain 
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accurate  numerical  computation  of  the  Fresnel  diffraction  integral.  Issues  associated  with 
sampling  are: 


1.  Discrete  samples  on  a  finite-sized  grid  requires  the  spatial  frequency  to  fulfill  the 
Nyquist-Shannon  sampling  theorem,  otherwise  aliasing,  rippling  and  smearing  will 
occur  in  the  spatial- frequency  domain  in  addition  to  virtual  periodic  replication  in  the 
spatial  domain. 

2.  Circular  shifting  of  the  samples  can  occur  in  the  spatial  domain  (i.e.,  samples  are 
time  shifted  so  that  the  first  sample  coincides  with  the  origin)  to  account  for  only 
positive  indices.  The  reordering  of  spatial  samples  means  that  the  samples  in  the 
spatial-frequency  domain  will  be  out  of  order,  too. 

3.  The  quadratic  phase  factor  inside  the  Fresnel  diffraction  integral  is  not  bandlimited 
and  hence  sufficient  spatial  resolution  is  required. 


In  addition,  the  geometry  due  to  laser  propagation  also  limits  the  observable  spatial  fre¬ 
quency  content  at  the  source  plane  when  using  Equation  5.3  for  single  step  propagation  due 
to  the  fixed  grid  spacing,  62  =  AAz/n6,,  at  the  target  plane. 

As  mentioned  earlier,  Equation  5.4  introduces  an  additional  scaling  parameter  by  manipu¬ 
lating  the  quadratic  phase  factor  which  is  inside  the  Fresnel  diffraction  integral  [17],  [18]: 


|r2  —  ri|2  =  r\  -  2r,  •  r2  +  r\ 


=  (r,  +  —  -  — 
\  “  m  m 


2ri  •  m  +  (fj  +  mr\  -  mr\ ) 


=  m 


m 


2  (  — )  •  ri  +r\ 
m  / 


+  1 


m 


r;  +  (1  -  m)r\ 


(5.7) 


r2  . 

—  -n 

m 


1  -  m 


m 


r;  +  (1  -  m)r\. 


Substituting  Equation  5.7  into  Equation  5.2  yields 
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E(r2)  =  ^  P  E(n)2*H(2HM‘?K+<'-"W](iri , 
J  — oo 


iAAz 

■iK  i(  k(l-m)  \  2 

js^lglkAz g  1  2Azm  )r2  f°° 
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r  1 
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E(ri)e 


'2AzH(m)  ri 


E(r2)  =  e 


l\  2A zm  )' 


xT~l  \T  E(n)e 


;fc(l-m)  |2 
2AZ 


|fl=ri 


Jri  , 


(5.8) 


Equation  5.8  has  the  exact  same  convolution  form  as  Equation  5.4  with  r2  replaced  by  r2/m 
and  A z  replaced  by  Az/m.  Since  the  quadratic  phase  factor  inside  the  Fresnel  diffraction 
integral  is  not  bandlimited,  the  waveform  that  is  to  be  Fourier  transformed  discretely  needs 
to  be  sufficiently  sampled  in  order  to  prevent  aliasing.  Hence  for  accurate  simulation  results, 
the  grid  spacing  and  the  number  of  grid  points  are  critical.  The  grid  spacing,  number  of  grid 
points,  and  propagation  distance  must  fulfil  the  following  constraints  based  on  the  Nyquist 
sampling  theorem  [12]. 


82  < 


AAz  -  D28 1 


D 1 


(5.9) 


Di  1  ft  |  AAz 
~  28x  +  282  +  28182  ’ 


(5.10) 


A  z\  AAz  /  A  z\  AAz 


(5-11) 


N  > 


AAzj 

min(did„)2  ’ 


(5.12) 


where  8 1  is  the  source  grid  spacing.  So  is  the  target  grid  spacing,  D\  is  the  source  dimension, 
Z>2  is  the  target  dimension,  and  A n  is  the  maximum  propagation  distance.  For  a  fixed  8 1, 
82,  and  Az,  the  number  of  grid  points  required  according  to  Equation  5.10  will  probably 
be  very  large  resulting  in  long  simulation  run  times.  In  order  to  shorten  simulation  run 
times,  multiple  partial  propagation  is  usually  preferred  [12].  Using  Equations  5.9  to  5.12, 
the  maximum  target  plane  spatial  resolution  is  ~  20  x  10_3m  requiring  at  least  512  grid 
points  when  a  one  meter  diameter  beam  director  with  source  plane  spatial  resolution  of 
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~  2  x  10_3m  is  focused  onto  a  target  five  kilometer  away  (assuming  that  the  maximum  spot 
size  is  0.5m).  The  maximum  propagating  distance  is  ~800m  and  hence  at  least  seven  partial 
propagations  are  required  to  avoid  aliasing.  Using  the  parameters  obtained,  the  beam’s 
transverse  diffraction  term,  in  Equation  4.5  for  a  spatially  combined  laser  can  be 

simulated  without  aliasing. 


5.2  Thermal  Blooming  Modeling 

The  additional  phase  incurred  due  to  thermal  blooming  described  by  Equation  4.5  will  be 
implemented  in  the  atmospheric  propagation  code.  The  change  in  temperature  obtained 
using  Equation  4.8  will  cause  a  change  in  the  ambient  refractive  index  according  to  Equa¬ 
tion  4.7,  which  alters  the  wavefront  as  shown  by  Equation  4.10  under  the  steady  state 
assumption  (i.e.,  aAT/at  =  0).  In  order  to  simulate  thermal  blooming  effects  during  laser 
propagation,  the  computation  process  is  iterative  because  the  change  in  phase  depends  on 
the  laser  intensity  that  is  absorbed  during  the  laser  propagation,  and  vice  versa  as  illustrated 
in  Figure  5.2.1. 


Figure  5.2.1:  Method  to  obtain  thermal  blooming  effects  for  laser  propa¬ 
gating  from  source  plane  to  target  plane  under  the  steady-state  =  0) 
assumption. 


Considering  absorption,  scattering,  and  thermal  blooming  effects,  Equation  5.8  becomes 


37 


The  steady  state  assumption  (i.e.,aAr/ar  =  0)  can  be  relaxed  to  incorporate  transients.  At  each 
time  step,  the  change  in  ambient  temperature  for  all  the  phase  screens  are  computed  using 
the  absorption  term  in  Equation  4.5  neglecting  the  conduction  effect  (assume  convection 
>>  conduction)  as 

[t+5ti  a 

AT}  =  /  —Idt  (5.14) 

Jt  pCp 

i- 1 

AT;  =  ^  AT*  +  AT}  (5.15) 

k= 0 

where  AT)  is  the  temperature  change  incurred  due  to  the  absorption  of  laser  power  during  the 
time  interval  from  t  to  t  +  8tj  and  A Tn  is  the  total  temperature  change  in  ambient  temperature 
at  time  tn.  During  the  same  time  interval,  the  convection  effect  due  to  wind  can  be  modeled 
by  shifting  all  the  phase  screens  according  to  the  wind  magnitude  and  direction.  This  whole 
process  is  repetitively  iterated  until  the  target’s  intensity  converges  as  illustrated  by  the 
Figure  5.2.2. 
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Figure  5.2.2:  Method  to  obtain  thermal  blooming  effects  for  laser  propagat¬ 
ing  from  source  plane  to  target  plane  under  relaxed  steady  state  condition. 


5.3  Incoherent  Beam  Combining  Modeling 

Equation  5.13,  which  models  coherent  sources  only,  can  be  modified  to  model  an  incoherent 
source  by  adding  a  random  phase  to  each  beamlet.  For  each  time  step,  the  random  phase 
to  each  beamlet  is  changed.  The  target  irradiance  pattern  is  then  obtained  by  averaging  the 
target’s  irradiance  pattern  over  several  iterations  until  it  converges. 


5.4  Genetic  Algorithm  Modeling 

The  arrangement  of  beamlets  at  the  beam  director  output  aperture  will  affect  the  amount 
of  laser  intensity  at  the  target  especially  for  coherent  beam  combination.  The  enclosure 
of  same  radii  and  intensity  beamlets  within  the  beam  director  output  aperture  is  similar  to 
a  circle  packing  algorithm,  which  is  a  very  complex  mathematical  problem.  Only  certain 
circle  packing  arrangements  have  been  proven  to  be  optimal  by  mathematicians  over  the 
past  few  decades.  Table  5.1  illustrates  circle  packing  arrangements  that  are  deemed  or  have 
been  proven  to  be  optimal.  From  Table  5.1,  14  or  more  beamlet  packaging  arrangements 
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are  conjectured  optimal  illustrating  the  complexity  of  circle  packing  algorithm.  For  a  given 
number  of  enclosed  beamlets,  the  ratio  of  the  beam  director  radius  to  the  beamlet  radius  is 
shown  in  Table  5.1. 


No  of  enclosed 

beamlets 

Ratio  of  beam 

director  radius  to 

beamlet  radius 

Diagram 

Remarks 

i 

i 

• 

Optimal  (Trivial) 

Fill  Factor:  100% 

2 

2 

$ 

Optimal  (Trivial) 

Fill  Factor:  50% 

3 

-2.154 

V 

Optimal  (Trivial) 

Fill  Factor:  64.66% 

4 

-2.414 

8 

Optimal  (Trivial) 

Fill  Factor:  68.64% 

s* 

-2.701 

* 

Optimal 

Fill  Factor:  68.54% 

6*  3 

‘Proved  optimal  by  Graham  in  1 968 

O 

Optimal 

Fill  Factor:  66.67% 

7 

3 

« 

Optimal  (Trivial) 

Fill  Factor:  77.78% 

8*1 

-3.304 

$ 

Optimal 

Fill  Factor:  73.28% 

g*i 

-3.613 

© 

Optimal 

Fill  Factor:  68.95% 

10*1 

* 1  Proved  optimal  by  Pirl  in  1969 

-3.813 

$ 

Optimal 

Fill  Factor:  68.78% 

ll*2 

* 2  Proved  optimal  by  Melissen  in 
1969 

-3.923 

$ 

Optimal 

Fill  Factor:  71.48% 

12*3 

-4.029 

« 

Optimal 

Fill  Factor:  73.92% 

13*3 

*  3  Proved  optimal  by  Fodor  in 

1969 

4.236 

$ 

Optimal 

Fill  Factor:  72.45% 

& 

14*1 

-3.304 

& 

Conjectured 

Optimal 

Fill  Factor:  74.74% 

Table  5.1:  Table  illustrates  how  to  pack  beamlets  of  same  radius  optimally 
at  the  beam  director  output  aperture.  Source:  [19]. 


Although  optimal  packing  of  same  radii  beamlets  is  important,  the  fill  factor  parameter  will 
also  affect  the  target’s  irradiance  magnitude.  Figure  5.4.1  and  5.4.2  illustrates  how  the  fill 
factor  can  be  different  for  a  six  beamlet  configuration:  one  with  six  equal  radius  beamlets 
that  is  deemed  to  be  optimally  configured,  and  another  that  improve  the  fill  factor  by  using 
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beamlets  with  two  different  radii.  The  fill  factors  for  Figures  5.4.1  and  5.4.2  are  66.67%  and 
79.18%  respectively.  Relaxing  the  "equal  radius"  constraint  may  allow  for  more  efficient 
packing  using  the  same  number  of  beamlets. 

In  addition,  a  genetic  algorithm  is  used  to  vary  the  intensity  of  each  beamlet  radii  group 
shown  in  Figure  5.4.3  to  investigate  if  the  target  intensity  will  increase.  A  genetic  algorithm 
is  a  probabilistic  optimization  method  which  involves  mutating  a  randomly  generated  pop¬ 
ulation  of  beamlet  intensities  towards  a  locally  or  globally  optimal  solution  after  multiple 
simulation  iterations.  Two  randomly  selected  intensity  sets  are  selected  and  propagated 
from  source  to  target  during  each  iteration,  then  the  set  that  produced  the  lower  target  peak 
irradiance  is  recombined  or  mutated. 
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Figure  5.4.1:  An  example  of  six  equal  radii  beamlets  enclosed  in  a  one  meter 
diameter  beam  director.  The  color  map  indicates  the  beamlet’s  intensity. 
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Figure  5.4.2:  Another  example  of  six  beamlet  configuration  comprising  of 
two  different  radii  beamlets  enclosed  in  a  one  meter  diameter  beam  director. 
The  color  map  indicates  the  beamlet’s  intensity. 
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Figure  5.4.3:  The  intensity  variation  of  different  radii  beamlet  configuration 
at  a  one  meter  diameter  beam  director  output  aperture.  The  color  map 
indicates  the  beamlet's  intensity. 
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5.5  Model  Validation 


The  best  method  to  validate  the  newly  developed  atmospheric  propagation  model  is  to 
compare  the  simulation  results  with  actual  experimental  data.  Due  to  time  and  resource 
constraints,  the  simulation  results  obtained  using  the  newly  developed  atmospheric  propa¬ 
gation  model  will  be  compared  with  the  results  obtained  from  WaveTrain  instead.  All  the 
WaveTrain  results  are  provided  by  Dr.  Pogue. 


WaveTrain  is  a  wave  optics  propagation  code  produced  by  MZA  Associates  Corporation 
available  to  U.S.  government  and  their  approved  contractors  only.  WaveTrain  simulates 
laser  propagation  by  solving  the  paraxial  wave  equation  4.5  for  coherent  sources  just  like 
the  newly  developed  propagation  model.  WaveTrain  accounts  for  thermal  blooming  effects 
in  a  self-consistent  transient  manner.  The  simulation  results  obtained  from  WaveTrain  are 
validated  with  experimental  data  and  hence  it  can  serve  as  a  benchmark  for  the  newly 
developed  atmospheric  propagation  model  [20].  The  beam  director  is  assumed  to  be  an 
off-axis  optical  system  on  a  stationary  platform  aiming  at  a  stationary  target  during  the 
entire  model  validation.  The  simulation  parameters  used  for  model  validation  are  listed  in 
Table  5.2. 


Parameter 

Value 

Laser  Wavelength 

3.8/rm 

Target  Distance 

5km 

Ambient  temperature 

300K 

Absorption  coefficient 

10-5m-1 

Scattering  coefficient 

10~5m-1 

Cross-wind  velocity 

lOm/s 

Laser  beamlet  type 

Uniform  distribution 

Table  5.2:  Parameters  used  during  model  validation. 
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5.5.1  Diffraction  Model 


The  diffraction  model  will  be  verified  using  both  analytical  means  and  WaveTrain.  The 
verification  is  performed  by  filling  the  entire  one  meter  beam  director  output  aperture  with 
a  single  lkW  laser  source  before  verifying  CBC  and  IBC  diffraction  effect.  Since  the  laser 
power  is  low,  the  laser  beam  propagation  path  should  not  be  affected  by  thermal  blooming. 
The  lkW  laser  source  at  the  beam  director  aperture  is  similar  to  a  point  source  illuminating 
a  one  meter  circular  aperture  and  hence  the  target’s  irradiance  pattern  should  comprise  of 
diffuse  circular  discs  known  as  Airy  patterns.  Since  the  target  distance  is  much  larger  than 
the  laser  beam’s  wavelength,  the  wavefront  arriving  at  the  target  can  be  approximated  by 
plane  wave.  The  target’s  two  dimensional  irradiance  pattern  for  a  one  meter  beam  director 
diameter  can  be  computed  using  the  Fraunhofer,  or  far-held  diffraction  solution  as 

/  =  /o(^7^)  ’  T  =  sin  0,  (5.16) 

where  k  is  the  spatial  angular  frequency  of  the  laser  beam,  J\  (y)  is  the  first  order  Bessel 
function,  and  I„  is  the  laser  intensity  at  the  beam  director  output  aperture.  The  J,<y>/y  term 
in  Equation  5.16  is  maximum  as  y  approaches  zero  (i.e.,  the  irradiance  is  maximum  at  the 
center  of  the  diffraction  pattern).  The  first  minimum  occurs  when  y  =  3.832  which  is  offset 
from  the  center  maximum  point  by 

Ar  x  Rx  6  =  Rx  1.22—,  (5.17) 

D 

where  R  is  the  target  distance  and  D  is  the  beam  director  diameter.  Based  on  Equation  5.17, 
the  offset  distance  from  the  central  maximum  to  the  first  minimum  should  be  2.32  cm  for 
a  target  distance  of  five  kilometer.  This  offset  distance  will  be  compared  to  the  simulation 
result.  In  addition,  the  normalized  intensity  cross  sectional  profile  will  also  be  compared 
with  analytical  ones.  Once  the  simulation  result  matches  the  analytical  computation,  the 
results  will  then  be  compared  with  those  from  WaveTrain.  After  this,  CBC  and  IBC 
diffraction  and  interference  patterns  will  then  be  compared  to  the  results  obtained  from 
WaveTrain.  A  bigger  target  spot  size  should  be  expected  for  IBC  as  the  spot  size  is  inversely 
proportional  to  the  beamlet  diameter  unlike  in  the  case  for  CBC  where  the  target  spot  size 
is  inversely  proportional  to  the  beam  director  diameter. 
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5.5.2  Thermal  Blooming  Model 


After  validating  the  diffraction  model,  the  entire  one  meter  beam  director  output  aperture  will 
be  filled  with  a  single  1MW  laser  source  to  verify  the  thermal  blooming  effect.  The  results 
obtained  will  be  verified  against  Wave  Train.  Upon  successful  thermal  blooming  validation, 
CBC  configuration  will  then  be  verified  against  WaveTrain.  If  the  IBC  diffraction  effect  and 
CBC  thermal  blooming  effect  matches  the  WaveTrain  results,  a  1MW  IBC  results  would  be 
hypothesized  to  be  right  since  WaveTrain  can  verified  for  a  coherent  source  only. 
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CHAPTER  6: 
Simulation  Results 


The  lethality  of  a  DE  weapon  depends  on  many  factors;  the  ability  to  melt  a  hole  on  the 
target  is  one  of  them.  This  capability  depends  on  the  target’s  material,  the  laser’s  far  field 
irradiance,  the  spot  size,  and  the  laser  dwell  time.  The  power  delivered  that  goes  into 
melting  the  target  area  is 

Pmelt  =  F  X  P bucket  ~  P/ossi  (6.1) 

where  Pbucket  is  the  laser  power  that  has  been  focused  into  the  so-called  "bucket"  that 
defines  target  area  (after  atmospheric  losses),  F  is  the  fraction  of  laser  power  absorbed  by 
the  target  (~0.2)  and  Pi„ss  is  the  power  loss  due  to  conduction  and  radiation  effects.  The 
power-in-the -bucket  parameter  is 


P  bucket 


I  dA, 


where  dA  is  an  area  element  on  the  target,  I  is  the  irradiance  on  the  target,  and  the  integral 
is  evaluated  over  the  bucket  area.  The  radiation  loss  is  approximately 


erMTLlt  ~  T 


melt  ambient ' 


where  s  is  the  blackbody  emissivity  of  the  target  surface,  cr  is  the  Stefan-Boltzmann  constant, 
A  is  the  surface  area  that  the  laser  beam  heats,  Tme[t  is  the  target’s  melting  temperature, 
and  Tambient  is  the  ambient  temperature.  The  conduction  loss  Pcon  shown  in  Figure  6.0.1  is 
approximately 


p  - 

1  con 


kArjnq(T,nelt  Tambient ) 


where  Aring  is  the  area  through  which  heat  is  conducted  away  from  the  target  area  (for  thin 
targets,  it  is  the  product  of  the  melted  hole  circumference  and  the  target’s  thickness),  k  is  the 
thermal  conductivity,  and  A.v  is  the  distance  of  the  induced  temperature  gradient  (~2cm) 
for  an  engagement  lasting  a  few  seconds. 
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Figure  6.0.1:  The  conduction  loss  directions  are  indicated  by  the  blue  arrows. 


The  melting  energy  required  is  material  dependent  and  is  given  by 


Qmelt  —  Q  heating  +  Q 


itransi 


(6.5) 


where  Q heating  =  Cpm(TmeU  -  Tambient)  is  the  heat  energy  required  to  reach  the  material’s 
melting  temperature,  Qtrans  -  mHf  is  the  energy  required  to  melt  the  material,  Cp  is  the 
material’s  heat  capacity,  Hf  is  the  material’s  latent  heat  of  fusion,  and  m  is  the  mass  of  the 
melted  target  volume. 


The  required  laser  dwell  time  is 

Tdwell 

The  average  irradiance  over  the  bucket  is 


Qmelt 

P melt 


T  _  Pbucket 

lave  ~  A,  ,  ’ 

n bucket 


(6.6) 


(6.7) 


where  Abucket  is  the  area  of  the  bucket.  The  minimum  irradiance  required  to  melt  a  certain 
volume  of  the  target  in  a  given  Tdweu  is  estimated  as 


P melt 
Abucket 
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(6.8) 


According  to  Equations  6.2  to  6.8,  the  DE  lethality  depends  highly  on  the  laser  dwell  time, 
the  minimum  target  irradiance,  and  the  power- in-the-bucket.  The  minimum  target  irradiance 
is  the  amount  of  power  per  unit  area  required  to  melt  the  target  while  the  power-in-the -bucket 
measures  how  much  laser  power  has  been  focused  onto  a  defined  spot  on  the  target. 

Consider  a  target  with  the  following  properties  in  Table  6.1.  In  this  case,  the  average 
irradiance  required  to  create  a  10cm  hole  on  the  target  within  eight  seconds  is  ~65MW/m2 
and  the  amount  of  power-in-the-bucket  required  is  ~0.5MW. 


Aluminum  properties 

Value 

Density, p 

2700kg/m3 

Thickness 

0.02m 

Melting  temperature,  Tmeit 

933K 

Latent  heat  of  fusion,  Hf 

400kJ/kg 

Heat  capacity,  Cp 

897J/kg-K 

Thermal  conductivity,  k 

237  W/m-K 

Ambient  temperature,  Tambient 

300K 

Emissivity,  s 

0.05 

Table  6.1:  Target  material  made  of  aluminum.  Source:  [21], 


Since  the  irradiance  and  power-in-the-bucket  parameters  are  crucial  for  DE  lethality,  they 
will  be  used  as  a  performance  indicator  during  validation  of  the  propagation  model.  Fur¬ 
thermore,  the  beam  director  fill  factor  affects  the  target’s  peak  irradiance  for  coherent  beam 
combination.  Different  fill  factor  configurations  listed  in  Table  6.2  will  be  simulated. 
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No.  of  enclosed  beamlets 

Remarks 

1 

single  beam 

3 

same  radii  beamlet  with  fill  factor  of  -65% 

6 

same  radii  beamlet  with  fill  factor  of  -67% 

6 

two  different  beamlet  radii  with  fill  factor  of  -79% 

13 

two  different  beamlet  radii  with  fill  factor  of  -84% 

Table  6.2:  Beamlet  configurations  used  for  simulation. 


6.1  Diffraction  Results 

Figure  6.1.1  shows  the  beam  director  and  target  irradiance  patterns  due  to  a  uniformly 
distributed  lkW  laser  source.  According  to  Figures  6.1.2  and  6. 1 .3,  the  peak  target  irradiance 
is  ~2MW/m2  and  the  first  null  distance  occurs  2.32cm  away  from  the  center  maximum.  The 
first  null  distance  obtained  from  simulation  matches  the  theoretical  analysis  performed  in 
Chapter  5.  Since  the  source  irradiance  pattern  is  symmetrical,  the  target’s  irradiance  pattern 
must  be  symmetrical  along  the  x-axis  and  y-axis  as  shown  in  Figures  6.1.2  and  6.1.3. 

Figure  6.1.1  also  shows  a  10cm  diameter  bucket  containing  -840W.  The  amount  of  power 
reaching  the  target  plane  is  -900 W.  This  result  is  consistent  with  Equation  4.1  where  -  1 00 W 
attenuation  of  the  laser  source  is  due  to  absorption  and  scattering  effects. 
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Figure  6.1.1:  The  irradiance  pattern  due  to  a  lkW  single  source  at  the  beam 
director  output  aperture  (left).  The  target  irradiance  pattern  with  a  10cm 
power-in-the-bucket  circled  in  white  (right).  The  color  map  represents  the 
intensity. 


x  10 


Y=0m 


x  (m) 


Figure  6.1.2:  The  cross-sectional  intensity  plot  at  the  target  plane  for 
Y=0cm  (left).  Magnified  view  of  the  cross-sectional  plot  showing  the  first 
null  point  (right). 
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Figure  6.1.3:  The  cross-sectional  intensity  plot  at  the  target  plane  for 
X=0cm  (left).  Magnified  view  of  the  cross-sectional  plot  showing  the  first 
null  point  (right). 


6.1.1  Coherent  Beam  Combination 

Figure  6.1.4  shows  the  beam  director  irradiance  pattern  created  by  three  uniformly  dis¬ 
tributed  coherent  beamlets  with  lkW  total  output  power.  According  to  Figure  6.1.5,  the 
target’s  irradiance  patterns  obtained  from  the  propagation  model  and  WaveTrain  are  essen¬ 
tially  identical.  From  the  cross-sectional  plot  in  Figure  6.1.7,  the  peak  target  irradiance 
is  ~1.3MW/m2  and  the  first  null  distance  along  the  x-axis  occurs  2.44cm  away  from  the 
center  maximum.  However,  the  first  null  distance  along  the  y-axis  occurs  2.69cm  away 
from  the  center  maximum  as  highlighted  in  Figure  6.1.6.  The  asymmetrical  source  irra¬ 
diance  pattern  causes  the  null  distance  to  be  different  along  the  different  axes.  Both  the 
propagation  model  and  WaveTrain  contain  -780W  within  a  10cm  diameter  bucket  (circled 
in  white  in  Figure  6. 1 .5).  The  amount  of  power  reaching  the  target  plane  is  ~900W  for  both 
the  propagation  model  and  WaveTrain. 
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Figure  6.1.4:  The  irradiance  pattern  due  to  three  coherent  beamlets  with 
lkW  output  power  at  the  beam  director  output  aperture. 
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Figure  6.1.5:  For  the  three  coherent  beamlets  with  lkW  output  power, 
the  target  irradiance  patterns  for  the  propagation  model  (left)  and  Wave- 
Train  (right).  The  10cm  diameter  power-in-the-bucket  is  circled  in  white. 
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Figure  6.1.6:  The  cross-sectional  intensity  plot  at  the  target  plane  for 
X=0cm  (left)  as  a  result  of  the  three  coherent  beamlets  case.  Magnified 
view  of  the  cross-sectional  plot  showing  the  first  null  point  (right). 


Figure  6.1.7:  The  cross-sectional  intensity  plot  at  the  target  plane  for 
Y=0cm  (left)  as  a  result  of  the  three  coherent  beamlets  case.  Magnified 
view  of  the  cross-sectional  plot  showing  the  first  null  point  (right). 


Figure  6.1.8  shows  another  beam  director  irradiance  pattern  created  by  13  coherent  beamlets 
of  two  different  radii  with  lkW  output  power.  According  to  Figures  6.1.9,  6.1.10,  and  6. 1.11, 
the  target  irradiance  patterns  obtained  from  both  the  propagation  model  and  WaveTrain 
have  nearly  identical  structures  and  peak  irradiance.  From  the  cross-sectional  plots  in 
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Figures  6.1.10  and  6.1.11,  the  peak  target  irradiance  is  ~1.7MW/m2  and  the  first  null 
distance  along  the  x  and  y-axes  occurs  ~2.32cm  away  from  the  center  maximum  which 
happens  to  be  the  same  as  the  single  source  configuration.  Both  the  propagation  model 
and  WaveTrain  contain  -770 W  within  a  10cm  diameter  bucket  (i.e.,  circled  in  white) 
in  Figure  6.1.9.  The  amount  of  power  reaching  the  target  plane  is  -900W  for  both  the 
propagation  model  and  WaveTrain. 


x(m) 

Figure  6.1.8:  The  irradiance  pattern  due  to  13  coherent  beamlets  of  two 
different  radii  with  lkW  output  power  at  the  beam  director  output  aperture. 
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Figure  6.1.9:  For  the  13  coherent  beamlets  case  with  lkW  output  power, 
the  target  irradiance  patterns  for  the  propagation  model  (left)  and  Wave- 
Train  (right). 
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Figure  6.1.10:  The  cross-sectional  intensity  plot  at  the  target  plane  for 
X=0cm  (left)  as  a  result  of  the  13  coherent  beamlets  case.  Magnified  view 
of  the  cross-sectional  plot  showing  the  first  null  point  (right). 


Figure  6.1.11:  The  cross-sectional  intensity  plot  at  the  target  plane  for 
Y=0cm  (left)  as  a  result  of  the  13  coherent  beamlets  case.  Magnified  view 
of  the  cross-sectional  plot  showing  the  first  null  point  (right). 


The  simulation  results  obtained  for  coherent  beam  combination  using  the  propagation  model 
are  summarized  in  Table  6.3. 
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Coherent  beamlet  configuration 

Power  in  a  10  cm 

diameter  bucket 

Remarks 

1 

~840W 

Symmetrical  first  null  distance  along  x  and 
y  axes  occurring  2.32cm  away  from  the 

center  maximum 

3  same  radii 

~780W 

Asymmetrical  first  null  distance  occurring 
2.44cm  and  2.69cm  away  from  the  center 
maximum  for  x  and  y  axis,  respectively. 

13  beamlets  with  two  different 

radii 

~770W 

Symmetrical  first  null  distance  along  x  and 
y  axes  occurring  2.2cm  away  from  the  cen¬ 
ter  maximum 

Table  6.3:  lkW  coherent  beamlet  configuration  results  obtained  using  the 


propagation  model. 


6.1.2  Incoherent  Beam  Combination 

Next,  the  following  incoherent  beamlet  configurations  will  be  validated  by  comparing  its 
peak  irradiance  and  spot  size  (i.e.  radius  measured  from  maximum  irradiance  to  ~14% 
of  peak  irradiance)  with  corresponding  results  for  the  coherent  configurations  using  the 
propagation  model. 

Figure  6.1.12  shows  the  irradiance  pattern  created  by  three  incoherent  beamlets  with  a  total 
output  power  of  lkW.  Comparing  with  Figure  6.1.5,  the  peak  target  irradiance  is  reduced  by 
a  factor  of  three.  This  reduction  factor  matches  the  number  of  enclosed  beamlets  according 
to  Chapter  3.  The  target  spot  size  shown  in  Figure  6.1.12  is  ~1.9  times  larger  than  the 
target  spot  size  shown  in  Figure  6.1.5.  Based  on  Table  5.1,  this  result  is  within  ~ 1 2% 
from  the  expected  ratio  (i.e.  2.15).  Although  the  peak  irradiance  is  reduced  by  a  factor  of 
three,  the  power-in-the-bucket  for  this  incoherent  configuration  is  only  ~20W  less  than  the 
corresponding  results  for  the  coherent  configurations. 
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Figure  6.1.12:  The  irradiance  pattern  due  to  3  incoherent  beamlets  with 
lkW  output  power  at  the  beam  director  output  aperture. 


Figure  6.1.13  shows  an  pattern  created  by  13  coherent  beamlets  of  two  different  radii  with 
lkW  total  output  power.  Comparing  with  Figure  6. 1 .9,  the  peak  target  irradiance  is  reduced 
by  a  factor  of  ~8.2.  Based  on  Chapter  3,  this  result  is  reasonable  as  the  reduction  factor 
should  be  larger  than  the  seven  equal  radius  beamlet  configuration.  The  target  spot  size 
shown  in  Figure  6.1.12  is  approximately  three  times  larger  than  the  target  spot  size  shown 
in  6.1.5.  The  power-in-the -bucket  for  this  incoherent  configuration  is  about  ~100W  less 
than  its  coherent  counterpart  due  to  lower  target  peak  irradiance  and  larger  target  spot  size. 


Figure  6.1.13:  The  irradiance  pattern  due  to  13  incoherent  beamlets  of  two 
different  radii  with  lkW  output  power  at  the  beam  director  output  aperture. 
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The  simulation  results  obtained  for  incoherent  beam  combination  using  the  propagation 
model  are  summarized  in  Table  6.4. 


Incoherent  beamlet  configura¬ 

Power  in  a  10  cm 

Spot  Diameter IBC 

IpeakCBC 

tion 

diameter  bucket 

Spot  Diameter CBC 

IpeakIBC 

3  same  radii 

-760W 

~1.9 

3 

13  beamlets  with  two  different 

-670W 

3 

-8.3 

radii 

Table  6.4:  lkW  incoherent  beamlet  configuration  results  obtained  using 
propagation  model. 


Based  on  both  the  analytical  computation  and  WaveTrain  results,  the  diffraction  effects 
simulated  in  the  propagation  model  are  accurate. 


6.2  Thermal  Blooming 

The  laser  output  power  at  the  beam  director  output  aperture  is  now  increased  to  1MW  so 
that  thermal  blooming  effects  become  significant.  Figure  6.2.1  shows  the  beam  director 
irradiance  pattern  created  by  a  uniformly  distributed  source  with  an  output  power  of  1MW. 
According  to  Figure  6.2.2,  the  target  irradiance  patterns  obtained  from  the  propagation 
model  and  WaveTrain  have  similar  structure  but  their  peak  target  irradiance  differs  by 
-28%.  This  difference  is  likely  due  to  different  thermal  blooming  methodology:  WaveTrain 
uses  a  transient  thermal  blooming  model  while  the  propagation  model  assumes  steady-state 
conditions.  Both  the  propagation  model  and  WaveTrain  irradiance  patterns  are  crescent 
shaped  and  bend  into  the  direction  of  the  transverse  wind  according  to  thermal  blooming 
theory.  Since  the  target  irradiance  patterns  from  the  propagation  model  and  WaveTrain  are 
shifted  away  from  the  origin  but  symmetrical  about  the  x  axis,  the  power- in-the -bucket  is 
computed  by  searching  for  the  maximum  power  that  is  within  the  10cm  diameter  bucket 
along  the  x-axis.  Figure  6.2.2  shows  a  10cm  bucket  containing  ~580kW  and  ~620kW 
computed  using  the  propagation  model  and  WaveTrain  results,  respectively.  Based  on  the 
power-in-the -bucket  values,  the  single  beam  configuration  can  melt  the  target  in  less  than 
eight  seconds.  Although  their  target  peak  irradiance  differs  by  -28%,  the  power-in-the 
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bucket  differs  by  only  ~7%. 


xIO5 


x(m) 

Figure  6.2.1:  The  irradiance  pattern  due  to  a  single  source  with  an  out¬ 
put  power  of  1MW  at  the  beam  director  output  aperture.  The  color  map 
indicates  the  intensity. 


Figure  6.2.2:  For  the  single  beam  case  with  1MW  output  power,  the  target 
irradiance  patterns  for  the  propagation  model  (left)  and  WaveTrain  (right). 


Due  to  time  constraints,  only  the  single  beam  configuration  will  be  simulated  and  anal¬ 
ysed  with  the  propagation  model  implementing  transient  thermal  blooming.  Comparing 
the  results  in  Figure  6.2.3  to  the  results  from  WaveTrain  in  Figure  6.2.5,  the  target  irra- 
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diance  pattern  are  almost  identical  in  structure  and  their  peak  target  irradiance  differs  by 
only  ~6°!o.  The  power- in-the-bucket  obtained  using  the  propagation  model  with  transient 
thermal  blooming  is  ~610kW.  This  results  only  differs  from  WaveTrain  by  ~2°7o.  From 
the  results  obtained,  it  can  be  concluded  that  the  propagation  model  with  transient  thermal 
blooming  produces  result  that  is  within  ~6°lc  of  results  from  WaveTrain.  The  runtime  for 
the  propagation  model  with  transient  thermal  blooming  is  roughly  similar  to  WaveTrain. 
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Figure  6.2.3:  The  target  irradiance  patterns  for  the  propagation  model  im¬ 
plementing  transient  thermal  blooming. 


Figure  6.2.4  shows  the  beam  director  irradiance  pattern  created  by  three  uniformly  dis¬ 
tributed  coherent  beamlets  with  1MW  output  power.  According  to  Figure  6.2.4,  the  target 
irradiance  pattern  obtained  from  the  propagation  model  and  WaveTrain  are  similar  in  struc¬ 
ture  and  their  peak  target  irradiance  differs  by  only  ~  1 5%,  which  is  much  less  than  the  single 
beam  configuration.  Figure  6.2.5  shows  a  10cm  diameter  bucket  containing  ~340kW  and 
~400kW  computed  using  the  propagation  model  and  WaveTrain  result,  respectively.  The 
power-in-the -bucket  for  this  beamlet  configuration  may  not  be  sufficient  to  melt  the  target 
in  eight  seconds.  This  beamlet  configuration  has  ~50%  less  peak  target  irradiance  and 
~40c/c  less  power-in-the -bucket  compared  to  the  single  beam  configuration.  A  single  beam 
that  has  equivalent  fill  area  for  the  three  beamlets  configuration  is  shown  in  Figure  6.2.6. 
According  to  Figure  6.2.6,  the  peak  irradiance  is  ~200MW/m2  and  the  power-in-the -bucket 
computed  is  ~450kW.  The  three  equal  radius  beamlets  configuration  has  —  1 8 c/c  less  peak 
target  irradiance  and  ~26%  less  power-in-the-bucket  compared  to  its  equivalent  fill  area 
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single  beam. 


Since  the  two  models  treat  thermal  blooming  differently,  it  is  not  surprising  that  they 
produce  different  results.  Nevertheless,  both  models  have  peak  irradiances  within  -30% 
and  power- in-the -bucket  within  ~15%  of  each  other.  The  performance  of  the  remaining 
beamlet  configurations  will  be  evaluated  using  the  propagation  model  except  when  the 
arrangement  of  beamlets  becomes  complex  (i.e.,  13-beamlet  configuration)  as  the  runtime 
is  significantly  faster  than  WaveTrain. 
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Figure  6.2.4:  The  irradiance  pattern  due  to  three  equal  radii  coherent  beam- 
lets  with  1MW  output  power  at  the  beam  director  output  aperture. 
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Figure  6.2.5:  For  the  three  equal  radii  coherent  beamlets  case  with  1MW 
output  power,  the  target  irradiance  patterns  for  the  propagation  model  (left) 
and  WaveTrain  (right). 


Figure  6.2.6:  A  single  beam  at  the  beam  director  aperture  with  the  same  fill 
area  and  output  power  as  the  three  equal  radii,  coherent  beamlets  configu¬ 
ration  (left).  The  target  irradiance  pattern  obtained  using  the  propagation 
model  (right). 


Figure  6.2.7  shows  the  beam  director  irradiance  pattern  created  by  six  uniformly  distributed 
beamlets  with  1MW  output  power.  According  to  Figure  6.2.7,  the  target  peak  irradiance 
is  ~150MW/m2  and  the  power-in-the-bucket  computed  is  360kW.  Based  on  the  power- in- 
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the-bucket,  this  beamlet  configuration  may  requires  more  than  eight  seconds  to  melt  the 
target.  This  beamlet  configuration  has  -56%  less  peak  target  irradiance  and  -38%  less 
power-in-the -bucket  compared  to  the  single  beam  configuration.  Although  the  fill  factor 
ratio  difference  between  this  case  and  the  three  equal  beamlets  configuration  is  only  -2%, 
their  relative  peak  irradiances  differ  by  -8%  and  the  powers-in-bucket  differ  by  -7%.  Thus, 
the  beamlet  arrangement  with  respect  to  the  transverse  wind  and  the  fill  factor  ratio  is 
important. 

A  single  beam  that  has  equivalent  fill  area  for  the  six  beamlets  configuration  is  shown 
in  Figure  6.2.8.  According  to  Figure  6.2.8,  the  peak  irradiance  is  ~210MW/m2  and  the 
power- in- the -bucket  is  ~460kW.  The  six  equal  radius  beamlets  configuration  has  -27%  less 
peak  target  irradiance  and  -22%  less  power-in-the -bucket  compared  to  its  equivalent  fill 
area  single  beam. 


x  (m)  X  (m) 

Figure  6.2.7:  The  irradiance  pattern  due  to  six  equal  radii  coherent  beam- 
lets  with  1MW  output  power  at  the  beam  director  output  aperture  (left). 
The  target  irradiance  pattern  with  a  10cm  power-in-the-bucket  circled  in 
white  (right). 
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Figure  6.2.8:  A  single  beam  at  the  beam  director  aperture  with  the  same 
fill  area  and  output  power  as  the  six  equal  radii,  coherent  beamlets  configu¬ 
ration  (left).  The  target  irradiance  pattern  obtained  using  the  propagation 
model  (right). 


To  further  investigate  the  effects  of  fill  factor,  another  six  beamlets  configuration  shown  in 
Figure  6.2.9  with  two  different  radii  and  a  total  output  power  of  1MW  is  tested.  According 
to  Figure  6.2.9,  the  target  peak  irradiance  is  ~170MW/m2  and  the  power-in-the -bucket 

computed  is  310kW.  Based  on  the  power-in-the -bucket,  this  beamlet  configuration  requires 
more  than  eight  seconds  to  melt  the  target.  This  beamlet  configuration  has  -51%  less 
peak  target  irradiance  and  -47%  less  power- in-the-bucket  compared  to  the  single  beam 
configuration.  This  beamlet  configuration  has  slightly  more  peak  irradiance  but  less  power- 
in-the-bucket  compared  to  the  six  equal  radii  coherent  beamlet.  The  decrease  in  power-in- 
the-bucket  may  be  due  to  the  asymmetrical  source  of  this  beamlet  configuration  with  respect 
to  the  wind. 

A  single  beam  that  has  equivalent  fill  area  for  the  six  beamlets  with  two  different  radius 
configuration  is  shown  in  Figure  6.2.10.  According  to  Figure  6.2.10,  the  peak  irradiance 
is  ~260MW/m2  and  the  power- in-the-bucket  computed  is  ~510kW.  The  six  different  radii 
coherent  beamlets  configuration  has  -39%  less  peak  target  irradiance  and  -39%  less  power- 
in-the-bucket  compared  to  its  equivalent  single  beam  configuration. 
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Figure  6.2.9:  The  irradiance  pattern  due  to  six  coherent,  different  radii  beam- 
lets  with  1MW  output  power  at  the  beam  director  output  aperture  (left). 
The  target  irradiance  pattern  with  a  10cm  power-in-the-bucket  circled  in 
white  (right). 


Figure  6.2.10:  A  single  beam  at  the  beam  director  aperture  with  the  same 
fill  area  and  output  power  as  the  three  beamlets  configuration  (left).  The 
target  irradiance  pattern  obtained  using  the  propagation  model  (right). 


Next,  13  beamlets  with  two  different  radius  and  a  output  power  of  1MW  shown  in  Fig¬ 
ure  6.2.11  is  tested.  This  beamlet  configuration  is  symmetrical  about  transverse  wind 
direction  and  has  fill  factor  slightly  more  than  six  beamlet  shown  in  Figure  6.2.9.  Accord- 
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ing  to  Figure  6.2. 1 1,  the  target  peak  irradiance  is  ~260MW/m2  and  the  power-in-the -bucket 
computed  is  410kW.  Based  on  the  power- in-the-bucket,  this  beamlet  configuration  may 
require  about  eight  seconds  to  melt  the  target.  This  beamlet  configuration  has  -25%  less 
peak  target  irradiance  and  -29%  less  power-in-the-bucket  compared  to  the  single  beam  con¬ 
figuration.  This  beamlet  configuration  has  the  closest  peak  irradiance,  power-in-the-bucket, 
and  irradiance  pattern  compared  to  the  single  beam  configuration  (i.e.,  best  peak  irradiance 
and  power-in-the-bucket  among  all  the  other  configurations).  This  result  illustrates  the 
importance  of  beamlet  symmetry  at  the  beam  director  output  aperture  and  fill  factor  ratio. 

A  single  beam  that  has  equivalent  fill  area  for  the  13  beamlets  with  two  different  radius 
configuration  is  shown  in  Figure  6.2.13.  According  to  Figure  6.2.13,  the  peak  irradiance 
is  ~280MW/m2  and  the  power-in-the-bucket  computed  is  ~520kW.  This  configuration 
has  -7%  less  peak  target  irradiance  and  -21%  less  power-in-the-bucket  compared  to  its 
equivalent  single  beam  configuration.  The  1 3  different  radii  coherent  beamlets  configuration 
also  has  the  best  peak  irradiance  and  power-in-the-bucket  among  all  the  other  configurations 
when  compared  to  their  equivalent  fill  area  single  beams. 

As  mentioned  previously,  this  beamlet  configuration’s  irradiance  pattern  will  be  compared 
to  WaveTrain.  Based  on  Figures  6.2.11  and  6.2.12,  the  irradiance  patterns  are  almost 
identical  and  hence  the  irradiance  pattern  for  any  beamlet  configurations  obtained  using  the 
atmospheric  should  be  roughly  identical  to  WaveTrain. 
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Figure  6.2.11:  The  irradiance  pattern  due  to  13  coherent,  different  radii 
beamlets  with  1MW  output  power  at  the  beam  director  output  aper¬ 
ture  (left).  The  target  irradiance  pattern  with  a  10cm  power-in-the-bucket 
circled  in  white  (right). 


Figure  6.2.12:  The  target  irradiance  pattern  with  a  10cm  power-in-the- 
bucket  circled  in  white  obtained  using  WaveTrain. 
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Figure  6.2.13:  A  single  beam  at  the  beam  director  aperture  with  the  same 
fill  area  and  output  power  as  the  three  beamlets  configuration  (left).  The 
target  irradiance  pattern  obtained  using  the  propagation  model  (right).  The 
color  map  indicates  the  intensity. 


6.3  Genetic  Algorithm 

The  genetic  algorithm  is  next  used  to  investigate  if  higher  peak  irradiance  can  be  obtained 
by  varying  the  intensities  of  beamlets  with  different  radii.  Therefore  the  genetic  algorithm  is 
configured  to  iterate  over  different  beamlet  intensities  to  optimize  the  peak  target  irradiance. 
The  six  coherent  beamlets  with  different  radii  shown  in  Figure  6.2.9  is  passed  to  the  genetic 
algorithm  to  search  for  maximum  peak  target  irradiance.  As  shown  in  Figure  6.3.1,  the 
peak  target  irradiance  is  ~190MW/m2  and  the  power-in-the -bucket  computed  is  340kW. 
Based  on  the  power-in-the-bucket,  this  beamlet  configuration  still  requires  more  than  eight 
seconds  to  melt  the  target.  This  beamlet  configuration  has  ~7%  less  peak  target  irradiance 
and  ~41%  less  power-in-the-bucket  compared  to  the  single  beam  configuration.  There  is  an 
improvement  of  ~6%  in  peak  irradiance  and  power-in-the-bucket  compared  to  the  constant 
intensity  beamlet  configuration.  Comparing  their  performance  with  their  equivalent  fill 
area  single  beam  configuration,  there  is  an  improvement  of  ~12%  in  peak  irradiance  and 
~6%  in  power-in-the-bucket. 
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Figure  6.3.1:  The  irradiance  pattern 
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Next,  the  genetic  algorithm  is  configured  to  optimize  the  peak  target  irradiance  for  the  13 
coherent  beamlets  with  different  radii  shown  in  Figure  6.2.11.  According  to  Figure  6.3.2, 
the  peak  target  irradiance  is  ~290MW/m2  and  the  power-in-the-bucket  computed  is  400kW. 
Based  on  the  power-in-the-bucket,  this  beamlet  configuration  still  requires  about  eight 
seconds  to  melt  the  target.  This  beamlet  configuration  has  ~  1 6%  less  peak  target  irradiance 
and  ~30%  less  power-in-the-bucket  compared  to  the  single  beam  configuration.  There  is 
an  improvement  of  ~9%  in  peak  irradiance  but  a  reduction  of  ~1%  in  power-in-the-bucket 
compared  to  the  constant  intensity  beamlet  configuration.  Comparing  their  performance 
with  their  equivalent  fill  area  single  beam  configuration,  there  is  an  improvement  of  ~1 1% 
in  peak  irradiance  and  reduction  of  ~2%  in  power-in-the-bucket. 
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Figure  6.3.2:  The  irradiance  pattern  due  to  13  coherent  beamlets  of  two 
different  radii  and  intensity  with  1MW  output  power  at  the  beam  director 
output  aperture. 
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CHAPTER  7: 
Conclusions 


A  new  steady-state  thermal  blooming  diffraction  code  has  been  developed  for  this  thesis. 
This  new  code  can  analyze  coherent  beam  combination,  providing  a  solution  that  produces 
higher  peak  target  irradiance  for  a  given  configuration  of  beamlets  with  different  radii  and 
irradiances,  and  computes  faster  than  WaveTrain  by  a  factor  of  three. 


7.1  Comparison  with  Analytical  and  WaveTrain  Results 

For  diffraction  effects,  the  propagation  model’s  results  concur  with  analytical  and  WaveTrain 
results  for  a  coherent  source  that  has  been  focused  onto  a  target  five  kilometer  away.  For 
thermal  blooming  effects,  the  peak  irradiances  obtained  from  the  propagation  model  is 
within  -30%  of  WaveTrain  and  the  power- in-bucket  computed  using  the  propagation  model 
is  within  ~15%  of  WaveTrain  too.  This  thermal  blooming  results  are  considered  to  be 
acceptable  for  these  extreme  cases  since  each  code  uses  a  different  thermal  blooming 
model. 


7.2  Beamlet  configuration  with  constant  irradiance 

Based  on  the  target  with  material  properties  in  Table  6.1,  all  beamlet  configurations  in 
Table  6.2  have  sufficient  irradiance  to  melt  the  target  but  insufficient  power-in-the  bucket 
to  melt  a  10cm  hole  within  eight  seconds  except  for  the  single-beam  and  13 -beamlet 
configuration.  In  fact,  the  13-beamlet  configuration  not  only  has  an  irradiance  pattern  that 
resembles  its  equivalent  fill  factor  area  single  beam  but  also  the  highest  peak  irradiance  and 
power- in- the -bucket  rating.  Therefore  the  13-beamlet  configuration  is  considered  the  best 
among  the  tested  configurations.  For  practical  implementation,  the  alignment  and  beam 
control  for  the  13-beamlet  configuration  will  be  complex. 
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7.3  Beamlet  configuration  with  varying  irradiance 


For  both  the  six-beamlet  and  13-beamlet  configurations,  there  was  a  slight  difference  in  the 
power- in-the -bucket  compared  to  the  corresponding  cases  with  constant  beamlet  irradiance. 
Using  the  genetic  algorithm  to  search  for  higher  peak  target  irradiance,  there  was  an 
improvement  of  ~10%  in  both  the  six-beamlet  and  13-beamlet  configurations.  In  fact,  the 
13-beamlet  configuration  has  ~8%  higher  target  peak  irradiance  than  its  equivalent  fill  area 
single  beam.  If  higher  peak  irradiance  is  required,  this  configuration  should  be  selected 
instead  of  the  same  intensity  configuration. 


7.4  Future  Work 

The  current  propagation  model  could  be  improved  to  incorporate  the  transient  model  de¬ 
scribed  in  Chapter  5  instead  of  using  the  steady-state  condition.  The  implementation  of 
the  transient  thermal  blooming  model  will  also  enable  the  current  propagation  model  to 
incorporate  turbulence  effects.  For  the  genetic  algorithm,  individual  beamlet’ s  intensity 
could  be  varied  asymmetrically  to  investigate  if  there  is  any  significant  improvement  in  the 
target  peak  irradiance  and  power- in-the -bucket  rating.  The  genetic  algorithm  should  be 
modified  to  search  for  the  highest  power-in-the -bucket  instead  the  peak  target  irradiance,  as 
that  is  a  more  significant  parameter  for  predicting  the  DE  weapon  effectiveness. 
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Appendix 


A.l 


Table  A- 1 .  Approximate  Laser  Power  Levels  Needed  to  Affect  Certain  Targets 


Multiple  perspectives  that  may  reflect  varyrjg  assumptions  about  BQ  and  ocher  factors 


Source 

Beam  power  measured  In  kMowaccs  (1<W)  or  megawatts  (MW) 

-10  kW 

Tens  of 
kW 

-100 

kW 

Hundreds  of  kW 

NW 

One  Navy 
briefing  (3010) 

UAVs 

SmelbocB 

Missiles  (sarong  at  SOO  kW) 

Another  Navy 
briefing  (3010) 

Short-range  operations 
against  UAVs.  BAH 
MANPADS  (SO  kW- 
lOOkW.  low  BQ) 

Extendedrange  operations 
against  UAVs.  BAM. 
MANPAD5,  ASCMs  flying  a 
crossing  path  (>100  kW. 
BQof-3) 

Operations  against 

- - 1-  *  ■  *  * 

supersonic,  mgnry 

maneuverable  ASCMs, 
vansonic  alr-eo-suriace 
missiles,  and  beUddc 
r*vssihes  (>l  MTW) 

Industry 
briefing  (3010) 

UAVs  and 
small  boats 
(SOkW) 

RAM  (IOO»  kW).  subsonic  ASCMs  (300 
kW),  manned  aircraft  (SOO  kW) 

Supersonic  ASCMs  and 
balls dc  missiles 

Defense 

Science  Board 
(DSB)  report 
(**») 

Swfece 

14  km 

Ground-based  air  and 
miisie  defense,  and 
cotmteripg  rockets, 
art  fiery,  and  mortars,  at  V 

10  km* 

"Baade  group  fefense  a'. 

S-20  kra  (l-J  MW) 

Northrop 
Grumman 
research  paper 

<MM> 

Soft  UAVs 
at  short 
range 

Aircraft 
and  cruise 

short 

range 

Soft 

UAVs  at 
long 
range 

Aircraft  and  cruise  missiles 
at  long  range,  and  artfiery 
rockets  (lower  hundreds  of 
kW) 

Artillery  she  Is  and  terminal 
defense  apinst  very  short 
range  beltetk.  missiles 
(tagher  lOOacf  kW) 

Table  A.l:  Approximate  Laser  Power  Level  Needed  to  affect  certain  targets. 
Source:  [1] 
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